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MÉTODOS VARIACIONAIS EM MECÂNICA DOS FLUIDOS 
E TRANSMISSÃO DO CALOR O 


RESUMO 


Numa curta introdução o autor refere as motivações que 
originaram a pesquisa de princípios extremantes em Fisica 
e o interesse da sua utilização sob o aspecto aplicado da 
solução analitica aproximada de problemas em Mecânica 
dos Fluidos e Transmissão do Calor. Analisa seguidamente 
através dum exemplo de condução do calor o Teorema de 
Prigogine para o minimo de produção de entropia, do qual 
mostra as limitações e a sua não unicidade, ral como do 
«potencial local» do mesmo autor e de Glansdor ff, como cri- 
tério de evolução em Física Macroscópia. Com o objectivo 
do cálculo aplicado é proposta uma formulação variacional- 
-parcial das equações do movimento e da conservação da 
energia para um fluído newtoniano em regime estacioná- 
rio. Às vantagens e inconvenientes dos métodos variacio- 
nais são apresentadas, 


o. INTRODUÇÃO 


por  J. DELGADO DOMINGOS 
Prof. 1,9, To 
N. E. E.M. do I,A.C, 


SYNOPSIS 


On a short introduction the motivations for the search 
of extremal principles in Physics are presented as well as 
their application to the aproximate analytical solution of 
problema in Fluid Mechanics and Heat Tranfer. Through 
an example in heat conduction the limitaticns of Prigoqne's 
Theorem on minimum entropy production are etressed. Its 
non unicity and that of the «local potential» of Glansdorff 
and Prigogine as a generul evolution criterium in Macros- 
copic Physics is referred. A new partial-variational for- 
mulation for the fluid motion 18 proposed and the energy 
equationa in steady state are presented, The advantages and 
disadvantages of variational methods for Fluid Mechanics 


Transfer problems are summarised, 


0.1 — A pesquisa de princípios extremantes originou-se de início em motivações teológicas e 


Maupertuis designou como «Princípio Metafísico» a formulação que apresentou para o mínimo de 
acção. A Euler — a quem se devem os primeiros trabalhos fundamentais sobre o cálculo das variações 
se atribui também a afirmação (!) de que «sendo a construção do Universo a mais perfeita possível, 
obra prima de um ser superior e infinitamente sábio, nada existe onde uma propriedade de máximo 
ou de mínimo se não revele. Por consequência, todos os efeitos podem ser deduzidos por métodos 
de máximos e mínimos tanto das suas causas finais como das mais eficientes». 

Sob uma ou outra forma, por motivos teológicos ou conceito de optimização nas leis da natu- 
reza, a pesquisa de princípios extremantes tem sido preocupação de físicos e filósofos desde sempre. 
Efectivamente, a formulação variacional em que figuram apenas um estado inicial e um estado final 


ocorre ooo 


(*) Duma lição proferida em 28 de Junho de 1967 no Imperial College of Science and Technology (Univer- 
sidade de Londres) a convite do Departamento de Engenharia Mecânica. 
(!) Em 1744, A citação é livre. 
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determinando o sistema a própria evolução através de um critério de eficiência máxima, de minimo 
de esforço, de optimização, em face de constrangimentos impostos, apresentou-se como fuga ao de- 
terminismo implícito nas equações diferenciais pelo qual toda a evolução fica determinada por con- 
dições iniciais e constrangimentos impostos. Sob este aspecto, a formulação variacional apresenta-se 
como mais fundamental e susceptível de conter o que equações diferenciais não poderão exprimir. 
Daí a tendência para lhes chamar princípios. 

Outro atractivo, e não menor, reside na capacidade de síntese, e a abertura para a exploração 
de domínios novos que a síntese por via de regra origina — a formulação Hamiltoniana da mecânica 
é sob este aspecto reveladora. 

Todavia, apesar do seu inequívoco interesse muitos dos «príncipios» não podem ser consi- 
derados como tais — em certos casos as condições de extremo não se verificam mas apenas as de 
estacionaridade; noutros, caso mais frequente, o extremo só existe relativamente a variações dum 
tipo particular (variações fracas) como é o caso do chamado «Princípio de Dirichelet», ou do «míni- 
mo de acção» na passagem duma configuração a outra [1]. 

Princípios extremantes só se conhecem, regra geral, para sistemas conservativos. Para sistemas 
puramente dissipativos em equilíbrio mecânico, cita-se frequentemente o «princípio do mínimo de 
produção de entropia», enunciado por Prigogine, o qual não pode todavia considerar-se como de 
validade geral dadas as hipóteses restritivas em que se baseia a sua dedução — trata-se efectivamente 
dum teorema no âmbito da Termodinâmica Irreversível Linear para o regime estacionário, que pos- 
teriormente referiremos em mais pormenor. 


0.2 — Para além das suas implicações de natureza física, a formulação variacional das equações 
fundamentais possui o maior interesse como meio de obter soluções analíticas exactas ou aproximadas 
de problemas cujo tratamento por métodos formais «clássicos» é totalmente impossível ou demasiado 
trabalhoso. É na perspectiva dos metodos analíticos aproximados que a presente contribuição deve in- 
serir-se. Deste modo, a própria natureza do problema físico em causa nos permite obviar a conside- 
rações de existência, continuidade, diferenciabilidade das soluções, tanto como à restrição a variações 
fracas porquanto a característica central deste tipo de tratamento reside na equivalência das equações 
de Euler à equação diferencial proposta, com condições aos limites que voluntáriamente se restrigem 
a valores impostos na fronteira. Muitas das restrições assim introduzidas podem levantar-se sem 
dificuldades de maior, apenas não as consideramos por agora porque a finalidade é realçar o inte- 
resse do método mais que o refinamento das suas possibilidades. 

Na sequência da exposição considera-se apenas o regime estacionário e o fluido incompressível. 
Por outro lado, embora se acentui que o objectivo central consiste no método analítico aproximado 
de solução das equações do movimento e energia para um fluido newtoniano, alguns considerandos 
serão introduzidos para mostrar a ligação com o teorema do mínimo de produção de entropia de 
Prigogine, e o do “potêncial local» de Glansdoff e Prigogine. 


1. A EQUAÇÃO DO CALOR E O TEOREMA DO MÍNIMO DE PRODUÇÃO DE ENTROPIA 


1.0 — O teorema do mínimo de produção de entropia demonstrado por Prigogine afirma que 
«para sistemas em equilíbrio mecânico, se forem válidas as relações de Onsager e os coeficientes 
cinéticos constantes, para condições na fronteira independentes do tempo, o sistema evolui para um 
regime estacionário que se caracteriza pelo mínimo de produção de entropia». 

Como exemplo, considere-se a condução do calor num sólido isótropo e não dilatável. A pro- 
dução de entropia por unidade de volume e de tempo é [3] 


q == grad (1) q 
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sendo q o vector densidade de fluxo de calor que satisfará a [4]: 


“T 
1) q= Lag grad (7) : 
pelo que 
1 2 
2) q == Ly grad (+) 


O coeficiente cinético Lag está ligado à condutibilidade térmica K por 
3) Laa = K T* 


Integrando 2) em todo o sistema termodinâmico obtém-se 


4) P= À Lag grad (1) AV 
“ V 


e P é a produção total de entropia. 
No âmbito do teorema do Prigogine Lay deve ser constante. Nessas condições, para constran- 
gimentos na fronteira independentes do tempo, P será mínimo em regime estacionário. 
Efectivamente, para condições na fronteira expressos por uma temperatura imposta e indepen- 
dente do tempo, a condição de estacionaridade de 4), exprime-se pela equação de Euler 


5) div (2 La grad a pai 


que por 2) se reduz a: 

(6) divqg=o0 

a qual constitui a equação do calor em regime estacionário e assume a forma habitual se usarmos 2) e 3) 
(7) div k grad T=0 


1.1. — Como se verifica, a obtenção de 7) foi artificiosa. O teorema exige constância de L«q, para 
a obtenção de 5) o que só é admissível para T=- constante em todo o sistema. Efectivamente, veri- 
fica-se experimentalmente ser K e não Laa que se mantém constante numa gama extensa de tempe- 
raturas. Por outro lado, admitindo Lga constante em todo o sistema introduz limitações importantes 
porque embora localmente o sistema (sub-sistema) se encontre pouco afastado do equilíbrio e seja 
portanto admissível a hipótese da termoestática quanto à aplicação num elemento do contínuo da 
sua equação fundamental (!), num sistema arbitráriamente extenso o afastamento global do equilíbrio é 
tão grande quanto se desejar. Nesta perspectiva, o teorema do mínimo de entropia é de mais restrita 
validade que a própria equação diferencial. 

1.2. — Posteriormente, Glansdorff e Prigogine [3], procurando a generalização do teorema referido 
para sistemas não lineares, foram levados a separar a produção de entropia devido à variação das forças 
generalizadas da proveniente da variação dos fluxos, concluindo que existiria o mínimo para a produção 
de entropia no que se refere à variação das forças. A produção correspondente aos dois efeitos seria 
igual em sistemas lineares com coeficientes cinéticos constantes, caso que corresponderia ao do teo- 
rema referido em 1.1. O desenvolvimento desta via conduziu os mesmos autores ao conceito de 


(') Conceito implícito em toda a formulação clássica das equações fundamentais dos sistemas contínuos. Na 
terminalogia considera-se a expressão de Gibbs em vez de equação fundamental, aliás equivalentes quanto às implica- 
ções da afirmação em causa, 
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«potencial local» considerado como critério geral de evolução em Física Macroscópica, isto é, os sis- 
temas evoluiram no sentido de minimizar o «potencial local» e o regime estacionário corresponderia 
ao seu mínimo. O conceito de potencial local foi recentemente interpretado à luz da teoria das 
flutuações por Prigogine num simpósio realizado na Universidade de Chicago em 1965 e dedi- 
cado a Termodinâmica, métodos variacionais e estabilidade. 

Não iremos discutir em pormenor estes trabalhos. Citam-se como exemplo de trabalhos recen- 
tes reveladores da importância do tema. A sua motivação foi de ordem física, mas também de natu- 
reza aplicada. A conclusão que se extrai é que o critério de evolução enunciado é artificioso e não único. 
Sob o aspecto das aplicações podem obter-se resultados equivalentes por via mais directa e realçar 
simultâneamente as «zonas* que carecem de esclarecimento sob o aspecto matemático ou físico. 


1.3. — Consideremos de novo a condução do calor num sólido que ocupa o domínio D com 
fronteira S, possuindo fontes de calor Q função dos parâmetros não especificados /;. 
A equação do calor será 


pe “= diva +Q 


PE a 
Multiplicando ambos os membros por — E -, integrando no domínio e usando o teorema de Green 
dt 


obtém-se : 


, 


. à be) - 
/ — EC e dV==| (q.n) e ds -J | osrad — + Q dV 
S » va 


Introduzindo a equação fenomenológica 


q=-—K grad T 


: Í 9 T ê 
e considerando na fronteira a temperatura constante = o ) obtém-se 
t f 


em fer (SO) av=< [| 5 ferad T)2—T o] iv fÊ ao 
o OR Ot 2 dt 


Resultando o sinal de desigualdade da condição de estabilidade intrínseca do sistema termodi- 


nâmico (c > 0) pois; > 0, o que confere ao primeiro membro um carácter intrinsecamente não 
positivo. 
Como se verifica, É varia monotônicamente com t pelo que o regime estacionário só é possível se 


O) ad O 


Verificando-se apenas o sinal de igualdade no caso de se atingir o equilíbrio termoestático em todo 
o sistema, 


dy 


| E pa JOB 
Como 4 > 0€e 573 < Oo regime estacionário o mao) so se atinge para é minimo, com- 
— é ) 


pativel com os constrangimentos. Portanto, o mínimo de 
(. .— K = 

8) =| — (grad T)2—-Q T |dV 
vi 2 ' 

é condição de estacionaridade do sistema. 


4 TECNICA N. 871 


Por outro lado 2 4==0, implica a Equação de Euler: 
9) div (K grad T) = Q 
que corresponde à forma diferencial da equação da condução do calor em regime estacionário. 


1.3.1. — Como se verifica, o mínimo de 6 e a equação 9) são perfeitamente equivalentes para 
os constrangimentos impostos. A dedução é directa e sem ambiguidades. Todavia & não corresponde 
à produção de entropia, nem a obtenção de 9) para o regime estacionário com fontes seria possível no 
ambito do teorema. 


K 
Em termos de produção de entropia, E representa uma fonte imposta (!) e E) (grad T)? a 


produção devido a condução do calor. 
Como $ > O implica que 


+ (grad T* => QT 


Dividindo ambos os membros por Tº obtém-se : 


K Q 
—— (gradT?) > — 
27º É = a 
o que exprime que a estacionaridade exige que a produção de entropia devida à condução exceda 
ou iguale o dobro da produção devida à fonte imposta. 


1.3.2. — Sob o aspecto aplicado, o interesse da formulação variacional reside na possibilidade 
de tratamento directo de 8) em vez da solução da equação diferencial 9). Para tal pode usar-se o método 
de Ritz,ou o de Kantarovich já exemplificado em trabalhos anteriores [6]. A sua convergência demons- 
tra-se para esse caso, e, embora a aplicação apenas exija a estacionaridade do integral, o facto de a 
solução exata corresponder ao mínimo do integral permite desde logo saber que a solução aproximada 
tende para a exacta por valores que lhe são superiores. 


2. O CONCEITO DE VARIAÇÃO PARCIAL 


2.0. — Consideramos até agora apenas o caso de K e Q não serem função do tempo. Todavia, 
ambos podem ser função de parâmetros, incluindo como caso particular a dependência do ponto . 6 
assume nesse caso a mesma forma e a obtenção de soluções exactas ou aproximadas processa-se de 
modo inteiramente análogo. 

Suponhamos porém que a dependência se verifica através da temperatura. Do ponto de vista 
físico, se é aceitável admitir que a variação de Q e Pé de segunda ordem relativamente a T, localmente, na 
vizinhança do regime estacionário, K e Q podem ser considerados constantes relativamente a peque- 
nas variações de temperatura. Deste modo, ambos podem ser considerados função do ponto se cal- 
culados em cada ponto em função da temperatura que nele se verifica no regime estacionário. 

Nestas condições, a formulação mantém-se e é susceptível de extensão a âmbito muito mais lato. 

A diferença fundamental quanto às restriçõs feitas em 1.0 e 1.1 relativamente a Lag reside em 
que neste caso a constância de Laq implicava a sua uniformidade em todo o domínio. No caso pre- 
sente essa uniformidade não é necessária porque pode existir dependência das coordenadas espaciais; 
a independência dos coeficientes fenomenológicos ou de Q apenas se refere a pequenas variações em 
torno do seu valor em regime estacionário. 


(') Proveniente da irreversabilidade associada a fluxos de outra natureza, por exemplo, carga eléctrica. Man- 
tém-se a terminologia habitual e não se consideram, ôbviamente, efeitos cruzados. 
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2.1. —Sob o aspecto matemático, o conceito anterior corresponde ao de variação parcial do 


integral. 
Retomando 1.3. e designando por To ==F(P) o campo das temperaturas no regime estacio- 


nário, pondo 
K =Ko [F (p)] OQ = Qo [A, F(P)] 
a formulação seria: 


Minimizar 


10) d - [ = (grad T*— Oo T| dV 
V 


e após variação (em à 6 = 0) introduzir a condição 
11) Ko=K (T) O, =0Q (hi, T) 


Como pode verificar-se, obter-se-á de novo a equação diferencial que corresponde ao regime esta- 
cionário. 


2.2. — Sob o aspecto aplicado, conhecer To significaria conhecer a solução do problema que se 
pretende resolver por tratamento directo do mínimo do integral. Todavia, duas vias são possíveis: 1) 
Solução iterativa — arbitrar To, calcular K==K (F (P)),e OQ = Q (P), minimizar o integral e verificar 
se o valor do To proposto corresponde ao que se obtém; se não corresponde recomeçar com uma 
nova função, eventualmente a determinada para T na primeira interacção — não existe porém demons- 
tração de que o processo convirja. 2) Solução directa: propor To e T com a mesma forma mas de- 
pendentes de parâmetros a determinar. Usando o método de Ritz tal consistirá em propor 


Substituindo T e T; em 10) resultará 


/ 


" [ K(Za vi) n ; on ) a | 
pre —— .grad | 2 Gy) — Oo (2 a A). «ci *f dV = 6 (a; ci) 
4 1 1 / 
Como K, e Qo não são variados, na minização de é os a; consideram-se como parâmetros, pelo 


que o mínimo de $ se exprimirá por 


d 


dci 


12) = 0 == f(a, c) e 


e a coincidência de Ty e T é assegurada exprimindo agora a condição 
13) ai = Ci 


com o que 12) e 13) constituem um sistema de 2n equações a 2n incógnitas o que permite a sua 


determinação. 
Os outros métodos directos do cálculo das variações podem usar-se de modo inteiramente 


análogo. 
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3. EQUAÇÃO DA CONSERVAÇÃO DA ENERGIA NUM FLUIDO NEWTONIANO 


3.0 — Na formulação de 10) considerou-se explicitamente a dependência de Q nos parâmetros 
dh, que podemos considerar as variáveis espaciais. Através do conceito da variação parcial foi pos- 
sível ainda consider Q função da variável dependente (T). Como a forma da função (P,T) não foi 
especificada, suponhamos que assume a seguinte forma: 


14) Q=pvóo—w. grad To 


se representarmos por W o vector velocidade, por 4 a função de dissipação : 


| 8 
ij o X1 ', Xj 
e u for a viscosidade do fluido suposto newtomiano e incompressível pode verificar-se que a estacio- 
naridade do integral 


15) a=[ [ [ Mo (gra DT us —w grad To) ) dV 


para valores de T impostos na fronteira tem como equação de Euler : 
16) div (K grad T) +- vó — w. grad To=0 


To corresponde ao campo das temperaturas e foi tratado como função do ponto em qw. Se após 
variação impusermos 


17) T=T 
resultará 
18) div (K grad T) + vó — w. grad T=0 


que é a equação da conservação da energia para um fluido newtoniano incompressível em regime 
estacionário. vw corresponde pois à formulação variacional parcial daquela equação. Poderá verifi- 
car-se ainda que a estacionaridade de 74 corresponde a um mínimo, se o regime estacionário for 
possível. 


3.1 — O maior interesse aplicado da formulação 15) reside na possibilidade do seu tratamento 
pelos métodos directos do cálculo das variações, tal como referido em 2.2. » possui todavia um sen- 
tido físico intrínseco que pode ser explorado, por exemplo, numa via paralela à expressa em 1.3.1. 


3.2 — No modo como foi introduzida, a expressão 14) aparece como uma hipótese à priori 
que se verificou satisfazer à formulação procurada. Todavia, a sua dedução pode fazer-se dum modo 
análogo ao considerado em 1.3, ou por via intuitiva partindo dos resultados anteriores. Efectivamente, 
considerando um elemento de volume, dois efeitos nos aparecem como equivalentes à «fonte de calor 
no sólido» — a dissipação viscosa, « é e um efeito contrário proveniente dum transporte de energia 
interna pelo fluido w. grad T. Esta forma provocaria uma dependência de T pelo que, tal como ante- 
riormente para Q, haverá que considerá-lo apenas função do ponto, e deste modo a forma final 
w. grad To. 


4 — EQUAÇÃO DO MOVIMENTO E DA CONSERVAÇÃO DA MASSA PARA UM FLUIDO 
NEWTONIANO INCOMPRESSÍVEL. 


4.0 — Por considerações análogas às que permitiram estabelecer as equações anteriores a formu- 
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lação variacional (parcial) que corresponde às equações de Navier - Stokes para um fluido incom- 
pressível em regime estacionário pode verificar-se corresponder á minimização do integral 


19) a= [[[ 
, j=1,2,3 


wº; e wº; não são sujeitos a variação; após variação do integral impor-se-á analogamente 


' | à quô 
ds a A (a Te) + Zo wi bd wº; |dV 
2 i d Xi Hj dxXj | 


20) Wº = wi 

Na expressão de 7: 
0, ON? o, á cá aço 

6=2(-— + -—)| — é a função de dissipação 

ij di 0X 

P = pressão hidrostática 

wi; == componente i da velocidade 

2; = força por unidade de massa 


4.1. Em 22, há quatro funções incógnitas e como tal sujeitas a variação: — wi, w3, W3,e P. É 
conhecido do cálculo das variações que a cada função corresponde uma equação Euler. Como pode 
verificar-se, a equação de Euler correspondente a w; corresponde à equação de Navier-Stokes para a 
direcção i. 

A equação de Euler correspondente a P reduz-se a 


isto é, à equação da conservação da massa para um fluido incompressível. 
Deste modo, na expressão variacional, 73, condensam-se as quatro equações do movimento 
para o regime estacionário. O caso particular de 723, em que se verifica o anulamento dos efeitos 


“ + . o Ww 
de inércia 


Xj 


Wi = o) corresponde ao chamado príncipio de Helmoltz [7] para os escoamentos 


uniformes. 


4.2. — Tal como referido para a equação da energia, o principal interesse de 7 reside na 
aplicação dos métodos directos do calculo das variações. 

No caso de 73, porém, se wi, wz e ws, forem propostos de modo que a equação de continui- 
dade seja satisfeita qualquer que seja o valor dos parâmetros ou funções arbitrárias, os termos em 
que figura P anulam-se idênticamente quando se faz a integração no domínio, o que reduz o trabalho 
do cálculo necessário. 

Procedendo do modo indicado obtém-se a função (aproximada ou exacta) para as velocidades 
as quais determinam P através das equações de Euler correspondentes a cada componente. 


5. INTERACÇÃO DA TRANSMISSÃO DO CALOR E DO MOVIMENTO 


5.0. —- Na forma como foi deduzido, 73 é independente da temperatura desde que | possa 
considerar-se constante. Deste modo, 72 determina o campo das velocidades que serão conhecidas 
para a obtenção de T através de 71. 

Sob O ponto de vista físico, 21e 73, exprimem realidades diferentes o que se nos afigura de 
importância — no primeiro caso é a irreversibilidade associada à transmissão do calor a caracteris- 
tica denominante. Por isso nele figura T como factor (de qualidade) na transferência da energia 
interna em causa. 73 tem apenas a característica duma energia, energia mecânica ou potêncial trans- 
formada irreversivelmente em energia interna. 21e 72, diferem dimensionalmente da variável T. Fisica- 
mente, a diferença pode justificar uma análise mais profunda por susceptível de esclarecer algumas 
das zonas obscuras na passagem da Termostática à Termodinâmica dos meios contínuos. 
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5.1. —Se as propriedades físicas são função da temperatura 721 e 73 ficarão simultâneos. Sob 
o aspecto da utilização dos métodos directos a via que se sugere é a consideração em 1, de w como 
parâmetro e em 72, de T. O sistema de equações que desse modo se obtém considerando os dois 
extremos em conjunto iguala o número de incógnitas, como se pode verificar, o que permitirá a sua 
determinação. 


6. REGIME TURBULENTO 


Toda a formulação anterior é aplicável ao regime turbulento na perspectiva «habitual» de o tratar 
a partir da introdução de um critério de média nas equações do movimento e da energia [8,...] e 
dos conceitos de «viscosidade turbulenta» e «condutibilidade térmica turbulenta» Sob o aspecto for- 
mal, uy e Kr serão função de ponto, cuja forma é conhecida à priori. Com ligeiras alterações, a de- 
pendência de x; nas derivadas da velocidade pode ser introduzida. 

A nosso conhecimento, métodos do tipo exposto nunca foram usados na consideração do 
regime turbulento. 


7. LACUNAS 


7.0. — A lacuna essencial da formulação apresentada reside nos aspectos formais que se ligam 
ao conceito de variação parcial. Da mesma dificiência enfermam as formulações Glansdorff e Prigo- 
gine para o «potencial local». Efectivamente, se o conceito pode ser aceitável sob o aspecto físico, 
sob o aspecto matemático a sua fundamentação é deficiente e as implicações de ordem prática cen- 
tram-se na estabilidade e na convergência das soluções analíticas aproximadas. 

Este é aliás um problema comum a muitos métodos numéricos onde a habituação parece ter 
feito esquecer que a lacuna existe e é importante. 

Presume-se que a formulação apresentada é susceptível de gerar pela consideração de formas 
apropriadas para as funções não variadas em 71 e 73, limites superiores e inferiores para a solução 
exacta. Efectivamente, com To (ou w;) escolhidos à priori obtém-se a formulação variacional exacta 
para um problema que não será geralmente o proposto mas que com ele se liga, e com o qual coin- 
cide quando o To escolhido igualar T. Quando To-=T a solução encontrada poderá ser um limite 
superior ou inferior. A consideração de ambas permitirá (se conhecidas) obter o erro máximo da 
solução aproximada. Para alguns casos particulares esta via já foi tratada com éxito o que parece 
validar a sugestão num âmbito mais geral. 


8. CONCLUSÕES 


8.0 A apreciação das potencialidades da formulação apresentada pressupõe certa familiarização 
com os métodos directos do cálculo das variações. Num intuito de sistematização resumem-se segui- 
damente o que pensamos ser, no momento actual, as suas principais vantagens e inconvenientes no 
que se refere a Mecânica dos Fluidos e Transmissão do Calor. 


8.1. Vantagens 


— Meio simples de obter soluções analíticas aproximadas e em muitos casos soluções exactas. 

— Estabilidade e convergência da solução, em geral. 

— Possibilidade de usar a experiência já existente noutros domínios (p. ex. Elasticidade) com 
estes métodos na solução de problemas em Mecânica dos Fluidos e Transmissão do Calor. 

— Métodos formalmente análogos para o tratamento de equações lineares ou não lineares. 

— Usando um computador, a obtenção de soluções requer usualmente menor tempo de cálculo. 
Além disso as solucões são analíticas e não apenas uma tabela de resultados. 

— Quando a experimentação física ou o conhecimento do fenómeno sugere o tipo de função a 
ser usado, o método pode dar uma solução muito aproximada com pouco labor de cálculo. 
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— Como o método fornece soluções analíticas, uma verificação pode ser feita em muitas cir- 
cunstâncias acerca da convergência e precisão do resultado pela consideração de valores 
extremos para os variáveis. 

— O método possui um significado físico intrínseco e como tal é susceptível de esclarecer o 
fenómeno em estudo ou facilitar a obtenção de soluções (p. ex. estabilidade de escoamentos, 
problema de Bénard...). 

— Pode usar-se como método puramente numérico a partir do cálculo numérico dos integrais 
que figuram na sua formulação. O esquema de aproximação é por vezes mais vantajoso que 
o que resulta do tratamento directo das equações diferenciais (esta via encontra-se práatica- 
mente inexplorada). 


s.2. Desvantagens 


— Sob um aspecto puramente formal as condições necessárias e suficientes de extremo não se 
encontram estabelecidas em termos gerais. 

— A obtenção de soluções exactas ou aproximadas requere em geral a consideração de muitas 
equações simultâneas (algébricas ou diferenciais ordinárias) a menos que se disponha pre- 
viamente de bons resultados experimentais ou do conhecimento do tipo de solução. 

O sistema de equações a considerar é usualmente não linear se as propriedades físicas do 
fluido não podem considerar-se independentes da temperatura. 


— Excepto num número restrito de casos não se conhece o grau de aproximação da solução 
obtida (tal como sucede nos métodos puramente numéricos). 
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CONDENSADOR DE SUPERFICIE 7.300 m2) PARA AGUA SALGADA, 
DO TURBO-GRUPO | DA CENTRAL DO CARREGADO 
peso total (c/ tubos de refrigeração ), cerca de 180 + 
comprimento = cerca de 14,2 m 
largura = cerco de 5,5 m 
altura = cerca de 7,2 m 
( Fabricado sob licença BROWN BOVERI CIE) 


PONTES ROLANTES, GUINDASTES E 
APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico 
TURBINAS HIDRÁULICAS — Fabrico segundo licença de A. C.M. de Vevey, S.A. 
TURBINAS A VAPOR -— —— Fabrico segundo licença de Brown Boveri, Cie. 
CALDEIRAS A VAPOR —————— — Projecto e fabrico segundo licença de Foster 
Wheeler, Co. 
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C. D. U. 624.391 


SISTEMA DIGITAL PARA TELEMEDIDA E REGISTO 
DE GRANDEZAS FÍSICAS (') 


RESUMO 


No presente artigo faz-se a descrição resumida dum sig- 
tema digital para transmissão ou registo de grandezas fi- 
sicas, apresentadas aos terminais sob forma analógica, Con- 
sideram-se também vs condicionalismos de concepção e de 
características finais, 


1 — INTRODUÇÃO 


Trata-se dum sistema digital que, por divisão 
no tempo, concentra 15 canais de informação 
numa única banda magnética ou numa só via 
telefónica. Compreende pois 2 terminais; um de 
emissão e outro de recepção. 

Todo o conjunto foi concebido no Laboratório 
Nacional de Engenharia Civil e destina-se à 
observação do comportamento de obras. 


2 — CARACTERISTICAS 


a) N.º de símbolos por palavra ......... 7 (128 
níveis) 

DD CORO assassina niaaa binário 
natural 


* CICLO MAIOR 


| 
| ,925ms CICLO MENOR 
| 


é ' 


por LUÍS VICENTE 


Engenheiro Especialista 
do Laboratório Nacional 
de Engenharia Civil 


SYNOPSIS 


The present paper is a quick description of a digita; 
system which is intended for transmission and recording 
of physical quantities, analoy in form at both terminal 
ends, Reference is also made how conception and final cha- 
racteristics were conditioned, 


v Tecnadoula passas paia atas resistên - 
cia — transistor (PNP) 
h) Codificação do tempo (dia, hora e minuto)... 
miRsSeSNDca soa NougTa por escantilhões de contacto 

i) Programação do registo magnético : 


— Arranque em instantes programados 
— Arranque por ocorrência de acontecimentos 
— Duração do registo 


A figura 1 resume as características da su- 
cessão codificada. 

Ocorre um impulso de peso nulo por cada 
ciclo menor, o qual evita a formação de falsos 
canais em branco. As sequências de 7 Os evi- 
tam-se nos canais de informação por limitação 
inferior das tensões de entrada. 


sims 


CANAL 15 | CANAL O | 
I SINCRONISMO 


Fig. 1 — Sucessão codificada 


cora CANAL 1 | CANAL 2 
c) Frequência de repetição ............... 4 kHz 
(8 kHz na leitura magnética) 
d) Campo de medida analógico ......... + 3V 
e) Banda de trequência por canal ...... O—3 Hz 
f) Sincronismo .......cecesesesescrereeneco «. Canal 


branco 


“o 


(*) A matéria do presente artigo foi apresentada 


de 1967 pela Standard Eléctrica. 
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3 — NATUREZA DA CONCEPÇÃO 


Procuraram-se harmonizar 3 qualidades que se 
consideram essenciais num sistema de memori- 
zação ou de telemedida de grandezas físicas: 
versatilidade, segurança e simplificação. 


ao 1.º Simpósio de Telecomunicações promovido em Abril 


un 


A maioria da instrumentação utilizada em En- 
genharia Civil traduz as grandezas a medir em 
tensões eléctricas. Considerou-se então que a 
tensão eléctrica seria a forma analógica que tor- 
naria o sistema mais versátil. Codificadores de 
posição são apenas utilizados na leitura do re- 
lógio-calendário, dispositivo que faz parte inte- 
grante do sistema. Mesmo assim, esta informação, 
traduzida em impulsos de tensão, é recodificada 
no conversor analógico-digital, único do terminal 
de emissão. 

A digitalização foi determinada pelo requisito 
segurança, o qual inclui estabilidade, precisão e, 
prâticamente, total eliminação de diafonia. Além 
disso, o método exclui a necessidade de cali- 
bragem. 

Admitiu-se que, em geral, as tensões de entrada 
têm valores suficientemente elevados para po- 
derem ser directamente codificadas sem prévia 
amplificação. É o que se passa com a maioria da 
instrumentação electrónica utilizada. Assim, não 
se consideram amplificadores, tanto na emissão 
como na recepção, evitando-se uma causa de 
erro e simplificando-se o sistema cujo propósito 
é apenas reproduzir na razão de 1:1. 

Sem se perder de vista o objectivo essencial, 
nem a natureza da informação a processar, as carac- 
terísticas finais do sistema foram em grande parte 
condicionadas a um critério de simplificação. Assim, 
por exemplo, tanto o número de canais (15), 
como o número de símbolos por palavra (7), per- 
tencendo à sucessão de termo geral 2"-1 permitem 
simplificar a obtenção das sequências de instru- 
ções. 

A sequência mais fácil, a partir de contadores 
binários, é de 2” instruções. Reservando-se uma 
para efeitos de separação, como sincronismo, por 
exemplo, restam sempre 2"-1 posições disponi- 
veis para enquadrar a informação. 

A largura de banda de 3 Hz por canal compa- 
tibiliza-se, pelo menos teoricamente, com uma 
frequência de amostragem de 6 Hz. Nestas con- 
dições, a frequência de repetição necessária seria 
aproximadamente de 770 Hz. No entanto optou-se 
por 4 kHz, a que corresponde uma frequência de 
amostragem por canal de cerca de 31 Hz. Desta 
forma, a mesma informação ocupa maior largura 
de banda. Mas isto permite a utilização de sim- 
ples filtros RC, tanto para entrada do terminal de 
emissão como para saída do terminal de re- 
cepção. 
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4 — TERMINAL DE EMISSÃO 


Limitar-nos-emos a uma descrição global do 
terminal de emissão pormenorizando apenas 
os aspectos que essencialmente particularizam o 
sistema. Será também essa a atitude adoptada 
quanto ao terminal de recepção. 


MONITOR 
LINHA 
Sd - erre 
E O q 
+ L | CONVERSOR 
hr A 
os > SEE: É 
a REGISTO 


(| MAGNÉTICO 


- | ANALÓGICO 
2a DIGITAL 


LEITURA 
DO PROGRAMAÇÃO 
TEMPO 


/ RELÓGIO 
ALE NDARI) 


4 / 
q / 
O ” 

Ed 


ne” 


Fig. 2 — Terminal de emissão 
Diagrama de blocos 


4.1 — Selector e conversor analógico-digital 


O gerador de sinais, a partir do qual são des- 
codificadas as diferentes instruções, tanto para o 
selector como para o conversor analógico-digital, 
é constituído por um multivibrador estável de 
16 kHz acoplado a um contador binário consti- 
tuído por 9 flip-flops. 

Os 2 primeiros flip-flops geram um ciclo de 2 
instruções (e duas pausas) que são utilizadas pelo 
conversor analógico-digital na determinação do 
valor de cado simbolo : 


1 — Impulso de ensaio 
2 — Impulso de decisão. 


Os 3 flip-flops seguintes geram 8 instruções 
que correspondem aos 7 símbolos de informa- 
ção por palavra, mais um simbolo separador e de 
peso nulo. 

Finalmente, os últimos 4 flip-flops condicionam 
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a sequência de 16 instruções destinadas ao cc- 
mando do selector electrónico. 

O canal O é utilizado para ensaio (monitor) do 
conversor analógico-digital, sem destruição da 
informação que eventualmente está a ser proces- 
sada através dos outros 15 canais. Na saída de 
registador (ou de linha) aquele canal é caracte- 
rizado, sempre, por uma sucessão de 7 Os. 

O comutador (selector própriamente dito) é 
constituído por 16 transistores (Fig. 3), cujos 
colectores são ligados às entradas e os emisso- 
res a um ponto comum: entrada do conversor 
analógico-digital. 


Q 
Co 
1 
C; 
|] 
] 
| 
| 
15 
Ca 


+6V 
Fig. 3 — Selector 


O estado dos transistores é, normalmente, de 
corte. O comando, que actua nas bases, impõe- 
-lhes sucessivamente o estado de saturação, ob- 
tendo-se, na saída comum, a tensão seleccionada: 
com uma diferença de poucos milivolts. 

A conversão analógico-digital processa-se por 
comparação da tensão de entrada com a tensão 
de saída dum conversor digital-analógico. O cir- 
cuito lógico que ensaia e decide vai condicio- 
nando este último à tensão de entrada e a partir 
dos símbolos mais significativos. Desta forma, 
tornam-se simultâneos os processos de quantifi- 
cação, codificação e transmissão. 

Os símbolos «marcados» no conversor digital- 
-analógico são sucessivamente memorizados em 
7 flip-flops. A tensão analógica é obtida no ponto 
comum duma sucessão geométrica (razão 2) de 
7 resistências que são electrônicamente comuta- 
das para as 2 tensões que constituem os limites 
de variação analógica:— 3 Ve + 3 V, 
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Cada um destes comutadores (elementos con- 
versores) têm 2 entradas que ligam directamente 
às 2 saídas de cada flip-flop. São constituídos 
por 4 transístores (Fig.4). Evita-se a saturação 
simultânea dos 2 transistores de comutação, pelo 
tempo de armazenamento daqueles que os ante- 
cedem e que, na condução, funcionam fortemente 
saturados. 


-3V 
= 5V - 
NES 4 
RX 
1 4 ; 
NX SAÍDA, 
y E” sd ANALÓGICA 
77/77) e 
+5V -5V 


+3V 
+5V 


Fig. 4 — Elemento conversor 


O sinal de saída do conversor analógico-digital 
é transformado numa sucessão de impulsos alter- 
nadamente positivos e negativos, eliminando-se 
assim a componente contínua e reduzindo-se 
práticamente a metade a largura de banda que 
é exigida ao meio de transmissão. 


42-— Codificação do tempo e programação 


Utiliza-se um relógio marca SAUTER, ZWE 
15 II, ao qual se adaptaram 3 escantilhões de 
contacto para codificação reflectida respectiva- 
mente do dia da semana (4 símbolos), da hora 
(6 símbolos) e do minuto (6 símbolos). 

A leitura do tempo é realizada por um anali- 
sador sequêncial a um ritmo de um símbolo por 
segundo. Um impulso negativo separa as dife- 
rentes palavras (Fig. 5). 


2" / 
.3v 
sv quit Ao fo, 
-3v | 
4 Dia 4 noRa £ o MINUTO 


Fig. 5 — Codificação do tempo 
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A programação realiza-se por ligação flexível 
das variáveis correspondentes aos diferentes sim- 
bolos (escovas de leitura e respectivos comple- 
mentos). O seu produto booleano, quando igual 
a 1, desencrava um temporizador CROUZET, ao 
qual se impõe o tempo de registo magnético. 
Esta operação pode ser igualmente desencadeada 
pela circunstância do valor duma de 4 grandezas 
físicas ultrapassar um dado limite. 


5 — TERMINAL DE RECEPÇÃO 


O terminal de recepção (Fig. 6) é essencial- 
mente constituido por: 

a) Regenerador — Rectificação de onda completa, 
para uniformização da polaridade dos impulsos. 
Regeneração de forma e de tempo. Obtenção 
duma sucessão uniforme para comando. 

b) Acumulador — Ordenação no espaço da infor- 
mação recebida. 

c) Conversores digital-analógicos — Conversão di- 
recta, a partir 
recebida. 

d) Selector — Geração de 15 impulsos de sincro- 
nismo, correspondentes aos 15 canais de infor- 
mação. Estes impulsos constituem instruções de 
leitura para os conversores digital-analógicos. 

e) Circuito de sincronismo de ciclo maior — Retorno 
a zero de selector por cada canal branco (7 Os 
consecutivos). 


do acumulador, da informação 


º 7) 


| 


| MAGNET a | E 4 


E RES RENT [ TRO | 
| > i ONv DIGITAL=-ANALÓGICO | | Ea e 


| o e PE O ar - 


| Sd (1 E i cê A pé 


pqp Ode 


ig. 6 — Terminal de recepção. Diagrama de blocos 


5.1 — Regenerador (Figs. 6 e 7) 


Como é necessária a obtenção duma sucessão 
de impulsos uniforme que se processe à fre- 
quência de repetição e em fase com os símbolos 
contidos no sinal recebido, o sistema deve incluir 
um oscilador sincronizado com esse sinal. A 
sucessão uniforme, por seu turno, contribui para 
completar o processo de regeneração. 
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ENTRADA 


—— DD —— S —b—» | 


E aa 7 


16 mm 

«uz [——S FF 2 | 
| 

Le DEI 


| 
| 
É q Lia d 
| SINAL 
| iba 4 VREGENERADO 
| 
Y 


(SUCESSÃON 
UNIFORME / 


Vig. 7 — Regenerador 


O oscilador é constituído por um multivibrador 
astável modulado em frequência. A tensão mo- 
dulante resulta da detecção de fase do sinal 
recebido, em relação ao sinal gerado. A cons- 
tante de tempo de detector absorve o «flutter» 
do registador magnético. 

Quando o sistema é utilizado em memorização, 
a velocidade de leitura é dupla da velocidade 
de gravação. A frequência de repetição é então 
de 8 kHz. 

Como a gravação se processa a 9,5 cm s !,o 
número de símbolos memorizados por centimetro 
de fita magnética é de 420, ao que correspon- 
dem, em média, 50 medidas realizadas com erro 
inferior a 1 “0. 

Na figura 8 mostram-se sucessivamente osci- 
logramas do sinal recebido, rectificado, parcial e 
totalmente regenerado. 


NSAER RANA 


1 AVATÁ 
Vw 


MMA ANIMA 


- - sm dd Aa . 


Fig. 8 — Regeneração. Oscilogramas 


5.2 — Conversores digital-analógicos 
Em cada um destes conversores, O sinal é me- 


morizado durante o ciclo maior. Assim, a tensão 
de saída conserva, durante esse periodo, o valor 
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que se observa no instante de leitura, obtendo-se 
uma resposta em «escada». Como o campo de 
variação analógica é igual ao do terminal de 
emissão, esta circunstância elimina a necessidade 
de amplificar, para se obter a reprodução à ra- 
zão de 1:1. 


6 — TECNOLOGIA 


O sistema é totalmente transistorizado. 

O circuito lógico básico é NEM (complemento 
da soma). As células de memória são constitui- 
das por flip-flops. 


A disposição mecânica dos circuitos impressos 
assemelha-se à construção tipo CS3 da Bell. Na 
figura 9 mostram-se os 2 terminais: de emissão 
e recepção. 


7 — CONCLUSÃO 


O sistema construído comporta-se de acordo 
com as características atrás mencionadas. Opera 
com segurança na gama de temperaturas am- 
bientes de O a 40º C. Constitui um poderoso 
auxiliar para a automatização de observação de 
obras de Engenharia Civil. 


Fig. 9 — Terminais de emissão e recepção 
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C.D.D. 624.344.5/94 


ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL (*) 


[— Breve nota mensal 


No conjunto do sistema as afluências do mês de 
Julho totalizaram 149 GWh, a que corresponde uma 
probabilidade de serem excedidas da ordem dos 77"/o 
e um coeficiente de produtividade de 0,73. 

As afluências acumuladas desde o início do ano 
hidrológico até ao fim de Julho têm uma probabilidade 
de serem excedidas da ordem dos 69" 


IH — Elementos gerais (GWh) 


| a) Mensais | Variação 
1966 | 1967 | Sh 
Produção hidráulica (Ph)... 423, 6 43º Ni + 4 
Produção térmica (P+)..... 0. DB EUT] — 
Produção total (PT). ..... 389,1 423,9, 1 8 
Energia recebida de empresas 
não pertencentes ao RNC (Er) 0,6 0,4 — 83 
Exportações (Ex)... .....'! 00 03 — 
Importações (1) ......... 0,0 1,2] — 
Saldo importador (Si). ..... V,O + 0,9 — 
Consumo em bombagem (Cp). | 0,0 0,0 O 
Produção para con- | (1) 
sumos perman. (Pep).... 351,8 3850 41 94 


Produção para con- 


sumos não perman. (Penp).. 12,7, T3A 0 

RR OtAR Ta Er Bo ESSE DR SUS 424,5! 458,14. 8 

Coeficiente de hidraulicidade o Oil 0,73 — 
NOTA 


(1) O aumento percentual da produção para consumos permanen- 
tes, tendo em conta a incidência dos domingos e dias especiais, é 
respectivamente de q9 e 10,8 €/g. 


II — Diagramas de carga dos dias característicos 


| 4.º feira; 


1966 1967 
Hredução hidraulica (Ph) MWh | 14969 16137 
Produção térmica (Pr) MWh' 0 O 
Produção total (PT) MWh | 14569 16137 
Trocas com [ Export. (Ex) MWh 0 () 
Espanha | Import. (1) MWh () 5 
“Cousum em bomb. hidroel. (Cb) Mun 0 o: 
Prod. para cons. perm. (Pop) MWh| 12527 13651 
Prod. para cons, não perm, (Penp)MWh| 2442 | 2491. 
TO PAL pr (Ex) win | 14969 16142 
mer máx. TIA s14 881 
se|p LEx) | Potência min. MW 402 435 
SE CTTOC Uriliz.da ponta horas 184 18,3 
E E | Factor de carga 0,77 0,76 
25 [Potência max, x. MW. E O 
Ge p Potência min. MW, 298 381 
Us “Pp Utiliz. da ponta horas 17,4 17,5 
| Fartor de carga | 0,73 0,73 
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IV — Energia armazunada nas principais albufeiras 


| No fim do mês 
Albufeiras: pigs pa 
wo | %() 


Alto Rabagão . «cc... o | 763,7 18,9 
Pitadela . o o uso o o oo] SEBR 69,1 
Venda Nova ; é o 1%». 0) 1098 80,3 
Salamonde . . +... E RI 26,1 J4,6 
Caniçada é E A Ens 29,9 90,3 
VIM” eso mos ms; 4 Sm sas 25,4 27,4 
E di mm Dn +] 263,0 71,5 
Castelo do Bode. . ., + ., 156,7 96,1 
Culhofrel sw cs Evo o 5,1 61,4 
Lagoa Comprida +. . . +... 27,6 (2) t,4 
Bqnta LUSA «pro rc esco 46,8 75,2 
PENNE » ds as sd ts . 10,6 82,2 
POvORO dé x GLS sé UE é 6,9 (3) 44,2 
Total com À, Rabagão +» |1623,6 76,8 

sem A, Rabagão. . .| 859,9 15,3 


NOTAS 


(1) Coeficiente de enchimento em energia das albufeiras. 


(2) Inclui 3,1 GWh armazenados em Vale do Rossim no início do. 
mês e 20 GWh no fim do mês. 


(3) Inclui 2,7 (Wh armazenados ne açude do Poio no início do mês 
e 1,7 GWh no fim do mês. 


(+) Elementos extraídos das estatisticas mensais do Repar- 
tidor Nacional de Cargas (R. N. C,) As produções e 
us consumos das empresas do R. N. OC, representam 
cerea de 940, dos totais do Pais. 
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O novo Nível automático 
de engenheiro 


WILD NAK 2 


com círculo horizontal 


jim, 


HEERBRUGG 


Para nivelamentos de média e grande precisão, 
com micrómetro de lâmina de faces plano- 
-paralelas para nivelamentos de grande precisão. 


Erro médio para 1 Km de 
nivelamento duplo: + 1,5 mm 


- Erro médio pera 1 Km de 
nivelamento duplo com micró- 
metro: + 0,4 mm 


Luneta de 30 >< de aumento 
com imagem direita de obser- 
vação 


Solicitem catálogo ou demonstração 
aos representantes exclusivos 


WILD PORTUGAL (DA. rc aas 
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INDUSTRIAIS TEXTEIS 


PARA OS VOSSOS MANCHONS, FLANELAS, TRANSPORTADORES, 
E TODOS OS TECIDOS TÉCNICOS 


EM FIBRAS NATURAIS OU SINTÉTICAS 


Consultem: 


PANAFEL: 


FÁBRICA NACIONAL DE FELTROS INDUSTRIAIS 
SOCIEDADE LIMITADA 


APART.: 9 ESCRITÓRIO E FÁBRICA 


TELEF.: 52093- 52094 OVAR 
TELEG: FELTROS dE Sud E ESTRADA DE S. JOÃO 


Empresa Wltramarina de Sondagens E Fundações, 1.º 


(SONDADORA) 


SONDAGENS geológicas e geotécnicas — CAPTAÇÕES de águas 
subterrâneas — REBAIXAMENTOS do nível freático — CONSOLIDA- 
* ÇÕES E ESTABILIZAÇÃO do solo — INJECÇÕES DE CIMENTO 
e outros produtos — INFRAESIRUTURAS de barragens, lúneis e pontes 
e FUNDAÇÕES de todos os lipos 


Única Empresa da especialidade com sede no Ultramar Português 


DELEGAÇÃO - LISBOA SEDE - LUANDA MOÇAMBIQUE 
R.S.Pedro de Alcântara, 1] Rua Serpa Pinto, 60 L.Marques - C.P.982 - Tel.741121 
Tel. 324693 e 213256 C.P.2178- Tel.4226 Beira - C.P. 1002 - Tel. 71497 
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C. D. U. 624.058.2: 539.374 
691.327 : 531,781,2 


PRINCÍPIOS TEÓRICOS DA DETERMINAÇÃO EXPERIMENTAL 
DIRECTA DE TENSÕES EM MEIOS VISCOELÁSTICOS. 
APLICABILIDADE AO BETÃO EM OBRA 


(Conclusão) 


6 — INTERPRETAÇÃO DOS RESULTADOS 


As expressões obtidas permitem o estabeleci- 
mento de novas bases teóricas para o projecto 
de tensómetros em meios que sofrem deforma- 
ções por fluência. 

Verifica-se que, num ensaio de fluência, as 
tensões suportadas por inclusões não completa- 
mente rígidas em meios viscoelásticos são com 
generalidade funções do tempo. Para que es- 
sas inclusões pudessem rigorosamente funcionar 
como tensómetros, deveria pois ter-se um factor 
de calibração (ou seja a relação entre a tensão 
a medir e a indicação fornecida pelo tensómetro) 
dependente do tempo. Na prática, poderá, porém, 
admitir-se a utilização de inclusões como tensó- 
metros, desde que a variação no tempo do factor 
de calibração se possa considerar desprezável, 
em face da precisão admissível para a medida. 
Através das expressões deduzidas na alínea an- 
terior, demonstra-se que, considerando como pa- 
râmetros o módulo de elasticidade e as propor- 
ções geométricas da inclusão, é teoricamente pos- 
sivel atribuir-lhes valores que tornem as variações 
do factor de calibração tão pequenas quanto se 
queira. 

Considerando o ensaio de fluência, o valor 
máximo do factor de calibração corresponde ao 
momento de aplicação da carga e o mínimo só 
será atingido ao fim de tempo infinito. Para uti- 
lizar uma inclusão como tensómetro, será neces- 
sário que a relação entre esses factores, que 
corresponde à relação entre as tensões supor- 
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tadas pela inclusão, não se afaste da unidade 
mais do que a percentagem de erro admissível. 

No que segue, analisaremos o quociente entre 
o máximo factor de calibração e o mínimo, ou 
seja, o quociente entre a tensão na inclusão ao 
fim de tempo infinito, e a tensão inicial. 

Vimos que, no caso de um meio viscoelástico 
regido pelas leis (3.1) e de parâmetros constantes, 
esse quociente era dado por (3.19) com o valor 
máximo (3.21). 

É fácil verificar que essas expressões tendem 
para a unidade se o valor de H/DE; tender para O. 
Isto é, na medida em que os problemas de sen- 
sibilidade da medida possam ser resolvidos, um 
tensómetro cilíndrico num meio com estas ca- 
racterísticas deverá ser uma inclusão extrema- 
mente rígida, ou uma inclusão de pequena rela- 
ção altura-diâmetro, ou ainda uma conjugação 
destes dois tipos, por forma a que o quociente 
H DE; não exceda um valor máximo conveniente, 
em face da precisão pretendida para a medida. 
Para cada grau de aproximação, aquele valor 
conveniente dependerá ainda, como se depreende 
da fórmula, dos valores de Em e de B, caracteris- 
ticas do meio indefinido. 


Se o fator tende para =, verificamos que 


a expressão (3.19) tende para 1 + Em B, o que 
será característico de um funcionamento como 
extensómetro. 

No caso de um tensómetro circular, a relação 
a considerar é (5.3) com o valor máximo (5.5), 
podendo pois extrair-se a mesma conclusão que 
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anteriormente em relação ao módulo de elastici- 
dade da inclusão para que possa funcionar como 
tensómetro. 

Com generalidade, poderíamos, para uma qual- 
quer forma de tensómetro, considerar a relação 


y + Em 
q =— E pi 
há J 1 = Em as | 
E; é + Em B 


em que > é um factor dependente da forma geo- 
métrica do tensômetro. 

No caso de um meio viscoelástico com pará- 
metros característicos dependentes da idade, e 
uma inclusão cilíndrica, a relação a considerar é 
intermédia entre a anterior, ou seja 


H Em EA nd » 
«D E; 2 
a E DX (39) 
«D E; 1+ Em Bo 2 
e 
H Em, 1-2? | 1+EmBo 
E DD Ei 2  1+ Em Bo 
% H ES Em 1 my2 nai 1+Em Bos 
aD E; 1 + Em Bs 2 1 + Ent Bo 
(5.19) 


desde que se tenha Boo < Bo. 

Demonstra-se facilmente que, dado que esta 
condição se verifique, o valor da expressão (5.19) 
é superior ao da expressão (3.19). Pode também 
notar-se que ambos os limites tendem para a 
unidade quando H/DE; tende para O, pelo que 
as conclusões anteriores são generalizáveis para 
este tipo de meio. No entanto, a velocidade com 
que o limite superior (5.19) tende para a unidade 
é agora dependente, não só de Em e Bo, como 
também de By. Se a relação (1+Eu Boo) /((1+ 
+ Em Bo) não é muito pequena, os valores das 
expressões (5.19) e (3.19) não são muito afastados. 
Caso contrário, o estudo feito deveria ser aper- 
feiçoado por forma a determinar rigorosamente 
a função da tensão no tensómetro ou aproximar 
entre si os limites inferior e superior a considerar. 

As expressões encontradas podem servir de 
base a uma crítica dos anteriores critérios de 
projecto de tensómetros em meios que sofrem 
deformações por fluência. 


No que respeita às bases teóricas originais da 
aplicação do método das inclusões nesses meios, 
verifica se que as expressões deduzidas lhes per- 
mitem dar uma justificação parcial. De facto, 
elas correspondem a assimilar o efeito das defor- 
mações de fluência a variações no tempo do mó- 
dulo de elasticidade de um meio elástico. 

Ora, num meio viscoelástico de parâmetros E e B 
invariáveis, sujeito a uma carga constante, atensão 
suportada por uma inclusão cilíndrica ao fim de 
tempo infinito é 


H a 
D " 2 í 
a(o) = — —— — 3 
Ha — Emo RR Semi * 
aD E 1+-EmnB 2 


(3.18) 


o que corresponde à tensão suportada por uma 
inclusão num meio elástico de módulo de elasti- 
cidade igual a Em /(1 + Em B), que é afinal o mó- 
dulo de elasticidade efectivo do meio viscoelás- 
tico ao fim de tempo infinito. Este resultado con- 
firma a interpretação teórica original do método 
das inclusões, neste caso particular, e legitima o 
critério para projecto de tensómetros, que consiste 
em determinar a relação E;/ Em a partir da qual 
a função (2.7) de E; / Em, ou seja, 


H + A (1 e 7) 
D 2 : 
G | =— EE e (G ums G 
H : Em. + a(1—v) — Qs 
D E 2 (6.2) 


se torne praticamente paralela à sua assíntota. 
No entanto, desde já se deve notar que a evolu- 
ção no tempo da tensão suportada pela inclusão 
cilíndrica no meio viscoelástico de parâmetros 
invariáveis não corresponde à sua evolução se 
inserida num meio elástico cujo módulo de elas- 
ticidade variasse no tempo de acordo com o mó- 
dulo de elasticidade efectivo do meio viscoelás- 
tico. Tal só se verificará no infinito, para o meio 
regido pelas leis (3.1). Por outro lado, é impor- 
tante notar que a tensão suportada pela inclusão 
num meio viscoelástico de parâmetros dependen- 
tes da idade (decréscimo com a idade da fluência 
específica) tende para um valor limite superior 
ao que corresponderia a um meio elástico de mó- 
dulo de elasticidade igual ao módulo de elastici- 
dade efectivo no infinito. Pode, finalmente, con- 
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cluir-se que, se a interpretação original do 
método das inclusões permitiu detectar os parà- 
metros de interesse no projecto dos tensómetros 
para meios viscoelásticos, e definir qualitativa- 
mente o seu sentido de actuação nos resultados 
obtidos, ela revela-se, por outro lado, insuficiente 
na definição das percentagens de erro consenti- 
das pelos tensómetros. As expressões apresenta- 
das neste trabalho são uma tentativa de progresso 
neste sentido. 

Como exemplificação, apresentam-se na fig. 2 
as curvas 9,9, q, eg, definidas respectiva- 


(9,94) 
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mente através de (3.19), (3.21), (5.19) e (6.2) em 
função da relação E;/Em para um meio viscoe- 
lástico em que se verifiquem 


1 + Em Bo == 4 
1+Em Bo 03 
1 + Em Bo 


e admitindo dois tipos geométricos definidos 
através de H/zx D = 0,2 e H/z D = 3. As relações 
escolhidas entre os parâmetros viscoelásticos fo- 
ram consideradas por se admitir que se possam 
verificar frequentemente no betão em casos típicos. 


q, do 


1,00 


0,80 


0,40 


0,00 
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Na fig. 3 estão representadas as mesmas cur- 
vas, mas agora em função de «x D/H, e mantendo 
constante o parâmetro E; / Em = 4. 

A função q, dá directamente, pela diferença 
entre a ordenada relativa a uma dada abcissa e a 
assintota da curva, ou seja a unidade, o máxi- 
mo erro consentido pelo tensômetro num ensaio 
de fluência se o meio é regido pelas leis (3.1). 
Por outro lado, as funções q, e q, enquadram 
o erro consentido pelo tensómetro num meio 
viscoelástico semelhante, mas em que o parâmetro 
B dependa da idade, decrescendo com esta. Veri- 
fica-se, a partir das curvas das figs. 2 e 3 que 
estes limites, em casos correntes, não serão muito 
distintos. Á incerteza no erro aumenta para va- 
lores baixos de E;'Em ou aD H. É possível que 
uma análise mais elaborada os permita aproximar 
se o interesse prático o exigir. A curva Qimáx é 


(9,94 ) 


1,40 


To NE l 
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bastante próxima de qr, para valores elevados de 
E;/Em e «D/H e é normalmente inferior a q',. É 
de notar ainda (fig. 2) que, para valo.es baixos 
de H/aD, os erros consentidos pelos tensómetros 
são muito pequenos, mas aumentam bastante 
quando a relação E;/Em é inferior ou pouco su- 
perior à unidade. Tal facto desaconselha o uso 
de tensómetros de compressibilidade semelhante 
a do meio, cujo projecto se fundamenta normal- 
mente na possibilidade de uma calibração directa 
do aparelho em máquina de ensaio e na diminuição 
das tensões secundárias de interacção com o 
meio. 

No que respeita à segunda interpretação do 
método das inclusões, baseada no conçeito de 
factor de dependência, as expressões agora dedu- 
zidas permitem concluir que aquele factor, no 
caso de meios viscoelásticos, é ainda uma função 
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do tempo. Para uma variação autogénia de volume 
do tipo da retracção, e que, portanto, se pode con- 
siderar como uma função do tempo, com derivada 
decrescente e de lei exponencial,o valor máximo 
(atingido ao fim de tempo infinito) da tensão 
secundária dada por (4.5) é relacionado com o 
valor fornecido para a extensão máxima pela teoria 
original (2.12) pelo factor de redução sy/ab, em 
que o significado das letras foi explicado na alií- 
nea 4. Pode concluir-se que a mesma relação 
existirá entre o valor máximo do verdadeiro factor 
de dependência para este caso eo valor original, 

É de assinalar o facto de o factor de depen- 
dência ser na realidade uma função das caracte- 
rísticas mecânicas do meio e da inclusão, e não 
somente geométricas, como na expressão original 
(2.12). Aquelas características são introduzidas 
através de a e b, por serem as raízes de uma 
equação cujos parâmetros são por elas condicio- 
nados. Nos casos que nos interessam (pensando 
no betão e no fenómeno de retracção) verifica-se 
que o factor de dependência diminuirá, embora 
não muito acentuadamente, com o aumento do 
módulo de elasticidade da inclusão [16]. 

Por outro lado, as deduções feitas permitem 
também concluir que o factor de dependência 
relativo às deformações autogénias de volume 
não pode ser considerado válido para as defor- 
mações de fluência. No entanto, se se pretender 
extrapolar o conceito de factor de dependência 
para o caso da fluência, dever-se-á admitir para 
ele uma nova expressão. Para extensão do con- 
ceito, o que se nos afigura mais lógico será con- 
siderar, para um ensaio de fluência, o quociente 
entre o acréscimo da tensão acusado pelo tensó- 
metro ao fim de tempo infinito e a tensão no 
momento de aplicação da carga. O factor de de- 
pendência para a fluência seria, portanto, dado 
por 
(6.3) 


Gi 


pelo que é fácil deduzir as suas expressões. No 
caso de meios viscoelásticos de parâmetros cons- 
tantes, teríamos 


H Em 1 — v* 
ec aes aD E; 2 | dest 1 
ai dm gt | 
«D E; 1+-Enm B 2 
(6.4) 
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Este factor de dependência corresponderia, por- 
tanto, ao máximo erro da indicação de um ten- 
sómetro cuja calibração fosse a do instante ini- 
cial. 

Faz-se notar que este factor de dependência 
tem uma nítida relação com o módulo de elasti- 
cidade do tensómetro, diminuindo com o aumento 
deste, como já se viu, e anulando-se para uma 
rigidez infinita. Este resultado invalida, portanto, 
a segunda interpretação do método das inclusões 
para a qual só a geometria da inclusão contava 
para efeitos de projecto de tensómetros. 

Para um meio em que o parâmetro B tendesse 
para infinito, teríamos o factor de dependência 
traduzido pela expressão muito simples 


2 H i Em 


——— O —uus 


ar(1—) D E 


(6.5) 


fmáx tn 


Devemos concluir, portanto, que, para O pro- 
jecto de tensómetros se poderá actuar, tanto na 
relação H/D como na relação Em/E;, redu- 
zindo-as a valores convenientes em face dos má- 
ximos erros pretendidos, para o que nos pode- 
remos socorrer das fórmulas apresentadas na 
alínea 5. 


7— COMPARAÇÃO COM OUTROS ESTUDOS 
SOBRE A INFLUÊNCIA DA VISCOELASTI- 
CIDADE NAS INDICAÇÕES FORNECIDAS 
PELOS TENSÓMETROS; INTERPRETAÇÃO 
E EXPLORAÇÃO DE RESULTADOS 


O estudo apresentado nas alíneas anteriores 
refere-se a meios viscoelásticos de características 
particulares. Designadamente, a lei reológica uti- 
lizada para a compressão simples 


G = Fm (t) o) 


E tt) = E Eis 


(7.1) 


traduz o comportamento de um corpo que, na 
literatura de língua inglesa, é frequentemente 
designado por «sólido standard» e de que o mo- 
delo reológico é constituído por um sólido de 
Hooke associado em série a um sólido de Kel- 
vin (ver fig. 4), pelo que também se designa por 
sólido de Hooke-Kelvin. 

A hipótese suplementar v == c.'º constitui uma 
outra opção, que nem sempre é adoptada no 
estudo teórico dos corpos viscoelásticos. De 
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facto, frequentemente admite-se que, sob carga 
constante, a componente volumétrica da defor- 
mação não varia com o tempo, isto é, aqueles 


Fig. 4 


corpos comportar-se-iam como elásticos para a 
deformabilidade volumétrica. Ter-se-ia, portanto, 
que 

Cj =3 K ii (7.2) 


em que o módulo de Bulk K é um valor cons- 
tante. 

No caso estudado, ao contrário, as originais 
das funções 


ú x 
pe a = 
2 Em" (1 + v) 
(7.3) 
K* A 


Re 3 Em* (1 — 24) 


permitem caracterizar o comportamento reoló- 
gico do corpo ao corte simples ou sob tensão 
hidrostática. 

A teoria desenvolvida na alínea 3 restringe-se, 
portanto, a um sólido «standard» com coeficiente 
de Poisson constante. Na alínea 5 os resultados 
obtidos são utilizados para enquadrar o que se 
passa num sólido do mesmo tipo, mas em que 
os valores das constantes características depen- 
dem da idade do material. Este caminho foi 
adoptado, uma vez que a preocupação funda- 
mental do trabalho era a determinação de ten- 
sões no betão, e se escolheram as expressões (5.3) 
fornecidas por Aroutiounian, como uma descrição 
analítica bastante ajustada do comportamento 
reológico daquele material. 
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As conclusões obtidas não podem portanto ser 
generalizáveis para todos os corpos viscoelásticos. 
Posteriormente à elaboração do nosso trabalho, 
tivemos conhecimento da existência de outros, 
em que a influência da viscoelasticidade sobre as 
indicações fornecidas pelos tensómetros é anali- 
sada. Um resumo destes trabalhos tem o maior 
interesse para efeito de complemento do estudo 
que se apresenta, e ainda para comparação de 
resultados. 

O primeiro é devido a Arlinger e Lundstrôm 
[21] (1961), e, conforme citação de Jan Hult |22] 
estuda um anel circular elástico embebido num 
corpo de Maxwell, e mostra que a evolução das 
tensões no tempo tem derivada decrescente de lei 
exponencial, e que a tensão do anel tende para um 
valor proporcional à tensão no corpo, mas inde- 
pendente das propriedades do meio envolvente. 

Jan Hult [22] (1963) apresenta uma tentativa 
de previsão do comportamento de tensómetros 
em diferentes tipos de meios viscoelásticos. 

Considera um modelo muito simplificado de 
distribuição da carga entre o meio e a inclusão 
em que estes funcionam sob a acção da carga S 
como os elementos de uma associação em para- 
lelo (ver fig. 5). 


| —+ inclusão 
m -— meio 


Fig. 5 


O autor classifica os tensômetros em passivos 
e activos (fig. 6'. Os tensómetros passivos podem 
ser esquemaáticamente representados por uma 
mola com uma certa deformabilidade. 

Sob a acção de uma carga CG, sofrem uma de- 
formação Ai proporcional à carga G através do 
inverso da sua rigidez c;, ou seja, 


Aj=—G (7.4) 
Ci 
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Se H é a altura do tensómetro e Q a sua secção 
transversal, é fácil verificar que a rigidez está 
relacionada com o módulo de elasticidade E; do 
tensómetro através de 


(7.5) 


a) tensómetro passivo 
b) tensómetro activo 
bi, b; — constantes dos tensómetros 


Fig. 6 


É de notar que, segundo este conceito, um ten. 
sómetro cilíndrico de elevada relação diâmetro- 
-altura será um tensómetro rígido, mesmo se o 
seu módulo de elasticidade for semelhante ao do 
meio envolvente. 

Os tensómetros activos têm outro princípio de 
funcionamento. A medida é feita fazendo-os 
actuar por uma pressão de óleo no interior tal 
que a sua deformação sob a acção da carga seja 
nula, Um tensómetro activo tem o mesmo efeito 
sobre o meio envolvente que um tensómetro 
passivo de rigidez infinita. Ambos os tipos de 
tensómetros podem portanto ser tidos em conta 
através do esquema da fig. 5, em que o modelo 
da inclusão é uma mola de rigidez c;. O resul- 
tado correspondente a um tensómetro activo é 
o que se obtém fazendo c; = cs. 

Quanto ao meio envolvente, o autor admite 
duas hipoteses: a viscoelasticidade linear e a vis- 
coelasticidade não linear. 
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No caso da viscoelasticidade linear, adopta o 
corpo de Burgers, constituído por uma associação 


em série de um líquido de Maxwell e um sólido 
de Kelvin (fig. 7). 


Fig 7 


Os elementos Hookeanos são caracterizados 
pela sua rigidez c e os dashpots pelas constantes r 
de proporcionalidade entre as forças aplicadas e 
as derivadas no tempo da deformação. 

A equação constitutiva deste corpo pode es- 
crever-se: 


À 4 8 À = ELS Crumrad 4 
Ta C E) ra 
+ Lig | 
LE T4 


O autor deduz que, por acção de uma força 
exterior constante S, a carga na inclusão, de 
valor Go quando introduzida no meio, tem uma 
evolução no tempo com derivada decrescente de 
lei exponencial, tendendo para o valor da pró- 
pria carga S, ou seja Goo =S (ver fig. 8). A carga 

GG 
2 


é atingida ao fim de um tempo dado por 
1 LN Tor 

ads pi A. ) Em; e E (7.6) 
Cc, Ci/ Tot Fa 
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o qual será mínimo se a rigidez da inclusão for 
infinita 


I ê Ta Ta) (7 7) 


O autor deduz ainda, utilizando a equação 
constitutiva geral com valores particulares dos 
parâmetros, as cargas na inclusão ao fim de tempo 
infinito (Goo) e os tempos 7 para os casos de 
um líquido não visceso (r, == 0), de um líquido 
viscoso (c,==r4=) de um corpo de Maxwell 
(r,==c0), de um sólido elástico (rs=r4 =), 
de um corpo de Kelvin (c, =r, = <<) e de um 
corpo de Hooke-Kelvin (r, = oc), 


s 


tos [ Go 
b) t=0* 


Fig. 8 


Finalmente, analisa os corpos viscoelásticos não 
lineares, admitindo que podem ser descritos muito 
aproximadamente pela equação constitutiva 


A=K.PMELK E (7.8) 


com n >m >1 


Considera dois casos de aplicação de força — a 
carga constante e a carga rapidamente variável — 
e conclui que, na maioria dos casós, é teoricamente 
possível a utilização como tensómetros de inclu- 
sões rígidas. 

Como o próprio autor faz notar, o modelo da 
fig. 4 é uma representação demasiado esquemá- 
tica de um corpo sujeito a carga e com uma in- 
clusão no interior, pelo que as conclusões extrai- 
das não podem ser consideradas válidas senão 
num sentido grosseiramente qualitativo; porém, 
o autor atribui-lhe a capacidade de deduzir as 
características mecânicas essenciais do sistema. 
O projecto definitivo de um aparelho concreto 
seria então necessáriamente completado por uma 
análise mais detalhada, em que fosse tida em conta 
a tridimensionalidade do sistema, a forma do 
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aparelho, etc., estudo esse do tipo do que é ex- 
posto no presente trabalho. 

Com base no estudo descrito de Hult, este autor 
e outros [23], (1966), elaboraram uma classifica- 
ção dos tensómetros em face das características 
do meio envolvente. 

Berry e Fairhurst [24] (1964) apresentam a 
resolução do problema de uma inclusão cilin- 
drica de comprimento indefinido inserida num 
meio viscoelástico, por uma forma que permite 
a previsão do comportamento do sistema num 
grande número de casos referentes às caracteris- 
ticas do meio. 

De facto, as deformações :; podem ser descritas 
a partir das tensões 7; através de funções do 
tempo Si; segundo a expressão geral 


“t 
e (t) = 5 (0) 2; (t) + | S'y (t=7) 2) (7) dr (7.9) 
“O 

Os índices i ou j, se iguais a 1, 2 ou 3, refe- 
rem-se a extensões ou tensões normais ; se iguais 
a 4, 5 ou 6, referem-se a distorsões ou tensões 
tangenciais. 

No caso elástico, as funções S; são indepen- 
dentes do tempo e podem escrever-se 


(=1,2,3) 


(i = 4,5, 6) 
(7.10) 


v 


a) (1, j=1,2,3) 
E J 


= 0 (,j=4,5,6) 


No caso viscoelástico, as funções S;; dependerão 

do tempo e podemos escrever (!) 

E;* = Si" Gj* (7.11) 

em que os símbolos * representam, como anterior- 
mente, a transformação de Carson. 

Os autores extraem conclusões de muito in- 
teresse com a única limitação de atribuir às fun- 
ções S11 e Si; as seguintes expressões, aliás bas- 
tante gerais 


(1) Os autores utilizam a transformação de Laplace ; 
porém, por simplicidade de exposição, continuaremos a 
utilizar a transformação de Carson. 
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de aços inoxidáveis, refractários e anti- 
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Culiesg d butA À “É Gli) (I-e—) dr 
ms (7.12) 
Ses (t) = as + bet + [o Go (7) . (I-e ) dr 


Sw e Si exprimem o comportamento do corpo 
à compressão simples ou ao corte simples; nas 
suas expressões, a; representa a componente ins- 
tantânea da extensão ou distorsão; b, a viscosi- 
dade; e Gi(7) são os espectros dos tempos de 
retardação da fluência. 

Retomemos, no entanto, o problema geral. Se- 
gundo o estudo de Savin [25] para o caso 
elástico, se Ni e N: são as tensões no meio 
actuando num plano normal ao eixo do cilindro, 
actuam neste tensões P; e P>, com as mesmas 
direcções de N; e Ns, e dadas por 


(Ni + Ny) 27 1T7 (N-No) 


nero a 
XT 


2 Ps 


f 


(7.13) 


em que o sinal * se refere à inclusão e y toma 
os valores 


3 - 47 
Pai É sed (7.14) 
é a 


conforme se trata, respectivamente, de um estado 
plano de deformação ou de tensão. 

Berry e Fairhurst interessam-se pelo caso da 
deformação plana e transformam a expressão em 


= a (Ni + N9)+b (Ni— Nº) (7.15) 
com 
sda == má 
— a+eu pr sEse 
u 


sendo c' e d' duas constantes da inclusão dadas 
por 


C=(u +? 
d' ;—1 (2, u -» constantes de Lamé) 
= | 


(7.17) 
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Nesta expressão, a solução está dada em termos 
das características do meio E e x que corres- 
pondem aos inversos de Si e Sis. 

A solução geral do problema no caso viscoelás- 
tico será, portanto, obtida através de 


gi = A* (Ny* + N$*) + B* (Ny* + N$*) (7.18) 


em que A* e B* são determinados a partir de aeb 
substituindo 1/E por Sw* e1/u por Ss”, obtendo-se 
facilmente 


Seo? (4Su* — Soo* 
25u* (Sc6* + c') 


DP: mm 


(7.19) 


ft —— 


| Seo* (4Su* — So6*) 
Sw? (7Su* — 2556*) + d'S u* 


Interessa, em particular, o caso em que Ni e N: 
são constantes. Teremos então Nj* = Ni; N4* = 
== Ns. Resultará, 


Pr(t) 


+ 3 
Ps (t) A (t). (Na + Na) * B(t) (Ni + N3) 


(7.20) 


em que A (t) e B (t) são as originais das funções 
A* e B*. 


Se agora adoptarmos para Su't) e Sw(t) as 
expressões atrás indicadas, virá 


Su =a + — + | — dY 
” go AH 
(7.21) 
S6* == de + a + | Se (O qu 
Pod Ep 


Neste caso, é fácil prever quais os limites para 
que tenderão, ao fim de tempo infinito, os valores 
das funções Pi (t) e Ps» (t'. Para isso, socorrem-se 
os autores do teorema de Tauber sobre as trans- 
formações de Laplace, que se pode escrever 


lim A (t) = lim A* (p) 


t— 00 p-—0 


(7.22) 


isto é, o limite no infinito de uma função do 
tempo é igual ao limite da sua transformada, 
função do parâmetro p, quando se faça tender o 
valor desse parâmetro para 0. 


No nosso caso, se temos bi, bs + 0, é fácil 
deduzir de (7.21) e (7.19) 


1 be 
A (00) =2— — — 
(ee) 2 b; 
(7.23) 
4 — b;/bi 
B (16) =——— 
7—2b;b:; 
Se by = bi == 0, 
oo: (Gi (2) 
Su *(p) = ay + —— dz 
u*(p) = as E pe 
o E fd (7.24) 
Soo *(p) = as + ni ii 
) » dE ps 


Se designarmos por Hi e Hs respectivamente 
os valores 


Hi = Su *(0)=a + N G (7) dt == Su (00) 
“a q 


e (7.25) 
Ho = Sw *(0) =as + | Go (7) dz = Sos (50) 


deduz-se fácilmente 


A (00) = At(0) = o (8Ht— Ho) 
2H: (Hs — c') 
H H (7.26) 
B (»0) ss B*(0) sm = 6 (4 à Hs) 


He (7 Hai — Ho) + d” Hi 


As conclusões a extrair das expressões (7.23) 
e (7.26) são curiosas. 

De (7.23) verifica-se que a tensão final na in- 
clusão elástica é independente das constantes 
elásticas da inclusão se o meio viscoelástico é tal 
que continua a fluir viscosamente tanto em tracção 
ou compressão uniaxial como em corte puro. As 
únicas constantes a saber são os coeficientes de 
que depende a deformação viscosa, ou seja by e 
bi. Contudo (exp. 7.26), se o meio viscoelástico 
não for deste tipo, a tensão final depende 
dos valores máximos das funções de fluência 
Su (00) =Hi ; Sos (o) = Hs; e, adicionalmente, das 
constantes elásticas c' = ()' + uy; ed =y i 
da inclusão. 

Os autores extraem deste estudo a conclusão 
geral importante de que a tensão numa inclusão 
cilíndrica pode variar apreciávelmente, mesmo se 
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não há variação do estado de tensão no meio 
envolvente; as variações de tensão determinadas 
por medições a longo prazo, utilizando inclusões 
elásticas devem, pois, ser consideradas suspeitas, 
a menos que tenham sido avaliadas as tensões 
devidas às deformações dependentes do tempo. 
O resultado expresso em (7.23) sugere ainda aos 
autores a observação de que deverá ser prati- 
cável, em rochas que exibam deformações viscosas 
apreciáveis, determinar as tensões N1 e N: no meio 
por observação dos valores limites P; e Ps para 
que tendem as tensões na inclusão. 

Estes resultados têm interesse para o nosso 
estudo, porque constituem elementos comple- 
mentares dos apresentados nas alíneas anteriores. 
De facto, alguns autores [26] admitem o modelo 
de Burgers para a descrição do comportamento 
reológico do betão. 

A função de fluência para a compressão simples 
deste corpo reológico pode escrever-se 


Ea tt) == Ê + At+B (Ie —!) 


m 


(7.27) 


Se admitirmos v= c.!* ou K =c.'s, este é um 
caso particular dos corpos tratados por Berry e 
Fairhurst, pois as funções Su (t)== Fm (t) e Sos (t) 
obedecem então às leis gerais (7.12). 

A solução (5.2) atrás apresentada refere-se a 
um problema semelhante ao tratado por Berry e 
Fairhurst, mas relativo a um estado plano de 
tensão, com Ni=N:==7, e definido num corpo 
de Hooke-Kelvin, com » == constante. Vejamos 
agora qual a solução se se tratasse de um corpo 
de Burgers ou de Maxwell, e admitindo ambas as 
hipóteses eventualmente adoptadas nas teorias 
da viscoelasticidade linear, ou seja v=-c'*, ou 
K = c!, 


No caso do corpo de Burgers, tem-se, de (7.27) 


1 A B; 
Ent ua =s JS E =çãs (7.28) 
Em P P +? 
A solução elástica 
2 E) 
EM 
(1 +ym) + (1) —S ço 
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transforma-se em 


* eum 2 E 
P it. — Vi 1 o 
E Masp ES 
Em 
(7.30) 
ui 2 Em 
é Med! 
vm) Ent — 
ou, por substituição, 
A B; 
E o pn 
pº = E p ie 
(1 + vm) TER 4 E | de 
m o P P+7 E; 
(7.31) 


O valor de p; (+=) será igual ao valor p; de 
(7.31) quando se fizer tender o valor do pará- 
metro p para zero, ou seja, como facilmente se 
deduz, 


Pp; (00) = perna (7.32) 
1 +m 

Este resultado é igualmente válido para um 
meio constituído por um líquido de Maxwell, 
cuja equação reológica é um caso particular de 
(7.27) para B = 0. 

Verifica-se que, similarmente à conclusão de 
Berry e Fairhurst, o valor final da tensão não de- 
pende das constantes elásticas da inclusão. É de 
notar que a expressão (7.32) é idêntica à expres- 
são (5.4) obtida para um sólido de Hooke-Kelvin 
quando o valor de B tende para infinito. Pela 
expressão (5.5) verifica-se que um tensómetro 
de rigidez infinita indicaria correctamente, em 
todo o momento, o valor exacto da tensão, após 
devidamente calibrado, visto que seria então 
pi (02) = pi (0). 

Se pretendermos admitir K == c!º, convém es- 
tabelecer a expressão (7.29) em função de E ek, 
para o que basta substituir o valor de vm pela 
sua expressão 


(7.33) 


resultando 


ES, 
| 
, 
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Pi=——————— a Õ 
3 1 —» 1 
— Fm + a 
2 E; 6 km 


Desenvolvendo F*m e determinando o limite 
para p= 0, deduz-se 


(7.35) 


4 
: (05) = —— 6 
Pi 3 


A tensão final na inclusão é, pois, proporcional 
à tensão aplicada, e não depende das constantes 
elásticas, quer da inclusão, quer do meio envol- 
vente. Este resultado coincide com o atrás citado 
de Arlinger e Lundstrôm. Verifica-se também que, 
neste caso, a tensão ao fim de tempo infinito 
numa inclusão rígida não é idêntica ao valor 
inicial da tensão, a não ser para o valor limite 
de vymm=0,5. Este é o valor que toma o coefi- 
ciente de deformabilidade transversal ao fim de 
tempo infinito sob carga. Então, o erro consen- 
tido per uma inclusão rígida em ensaio de fluên- 
cia em compressão num corpo de Burgers com 
módulo de Bulk invariante seria 


1 + (vmo em 1 + (Ym)o 
1,5 1 + (vm)oo 


ou seja, a relação entre a tensão supor rtada pelo 
tensómetro ao fim de tempo infinito e a tensão 
no momento de aplicação da carga, em que (mo 
e (vm). seriam os valores dos coeficientes de 
Poisson instantâneo e do coeficiente de deforma- 
bilidade transversal ao fim de tempo infinito. 

Conclusões semelhantes se podiam tirar rela- 
tivamente ao problema da inclusão cilíndrica no 
meio indefinido. 

Assim, se, por exemplo, este meio tem um 
comportamento em compressão simples regido 
pela equação reológica (7.27), sendo, além disso, 
vy==!€, obteriamos 


o (=0| 5 + 20-82 | 
P B7YP 
E de fue 
= Em p+y 
—p + ra(i=v)= (P ra SA 
DE; m p+y7 
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resultando, da aplicação do teorema de Tauber 
que 


Pode verificar-se facilmente de (2.7) que esta 
tensão equivale à suportada no instante t==0 
por uma inclusão completamente rígida, pelo que 
este tipo de inclusão poderá funcionar neste caso 
como tensómetro. 

Se a equação reológica em compressão simples 
é ainda dada por (7.27) mas, agora, é invariante 
o módulo de Bulk, obteriamos, após substituição 


em (2.7) de » por (> og ). e de — por Em”, 


= (5 + cas) Ea de ei Rs 
12 K 72Kº? 


3 . PRB 1 af 
——— (X E *2 + = 4- =) Em” e 
8 á ! DE 12kK/ 72K? 


e, por aplicação do teorema de Tauber, após subs- 
tituição de F*m (p) pela sua expressão (7.28), vem 
s Ha 

3 D ax 


Er (00) == 


+1 


tensão independente das características elásticas, 
quer da inclusão, quer do meio envolvente, e que 
corresponde à tensão suportada no instante t==0 
por uma inclusão cilíndrica completamente rígida 
se o meio envolvente tivesse um coeficiente de 
Poisson vm ==0,5. Nos outros casos, a inclusão 
rígida consentiria um erro dado por 


5 H 1 


1 
ET (00) o 3 D x E 
7 (0) a i. “+17 
D z(1—vmno) 
H 1 
2 — —+1 
D a«(1—-Ymco)r 
H 1 
Z2— + ———— - 1 
D «(1 —vmo) = + 


Este erro será mínimo se H/D tende para zero. 
Uma inclusão completamente rígida e com 
uma elevada relação diâmetro-altura seria, por- 
tanto, um bom tensómetro num meio deste tipo. 


28 


Os resultados obtidos para estes exemplos pa- 
recem permitir generalizar a conclusão de que 
inclusões completamente rígidas inseridas em 
corpos com comportamento de Maxwell ou de 
Burgers em compressão simples e com coeficiente 
de Poisson invariante suportam permanentemente 
a mesma tensão proporcional à tensão constante 
aplicada ao meio envolvente; porém, se é o mó- 
dulo de Bulk que se mantém invariante, a tensão 
suportada pela inclusão varia desde o seu valor 
instantâneo a um valor final que corresponderia 
ao valor instantâneo se o coeficiente de Poisson 
do meio envolvente fosse igual a 0,5. No caso 
de inclusões cilíndricas, será possível actuar sobre 
a relação altura-diâmetro, por forma a que o 
desvio entre o valor final e o instantâneo seja 
tão pequeno quanto se queira, o que se consegue 
baixando o valor daquela relação. 


8 — APLICABILIDADE AO BETÃO EM OBRA. 
CONCLUSÃO 


Nas alíneas anteriores, houve oportunidade 
de estudar o comportamento de inclusões elás- 
ticas em meios viscoelásticos lineares de dife- 
rentes tipos sujeitos a forças exteriores mantidas 
constantes. Designadamente, foram considerados 
o sólido standard ou de Hooke-Kelvin, o corpo 
de Burgers, o corpo de Mawell e o sólido stan- 
dard em que os parâmetros característicos depen- 
dessem da idade do material. Verificou-se que, 
nos casos em que o coeficiente de deformalidade 
transvsersal do meio não dependa do tempo de 
actuação de carga, as inclusões completamente 
rígidas podem teóricamente ser utilizadas como 
tensómetros precisos, o que constitui a base de 
uma teoria de determinação experimental directa 
de tensões em meios visco-elásticos. 

A conclusão exposta nem sempre será exten- 
siva aos casos em que aquele coeficiente não seja 
invariante; porém, no caso normalmente consi- 
derado de o módulo de Bulk ser invariante, é 
previsível o máximo erro consentido por um 
tensômetro completamente rígido, o qual não 
será grande se o coeficiente de Poisson do meio 
não for excessivamente baixo (por exemplo, 
inferior a 0,2). Foram ainda referidos estudos 
em que, embora adoptando-se esquemas muito 
simplificados de comportamento estrutural do 
conjunto meio-inclusão, se conclui pela validade 
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de utilização de tensómetros rígidos em meios 
viscoelásticos não lineares bastante gerais [22]. 

O betão é um material de comportamento 
reológico muito complexo, que se não pode filiar 
rigorosamente nos esquemas simples de visco- 
elasticidade linear acima citados. No entanto, a 
generalidade com que actualmente é posta em 
teoria a viabilidade de determinação experimental 
de tensões em meios viscoelásticos, parece legi- 
timar a sua aplicação ao betão. Por outro lado, 
são correntemente aceites como bastante apro- 
ximadas para fins de aplicação expressões de 
comportamento reológico deste material corres- 
pondentes ao sólido standard de parâmetros 
variáveis com a idade, e coeficiente de deforma- 
bilidade transversal invariante [17]. A análise 
conduzida na alínea 5, mostra que as conclusões 
gerais atrás citadas são ainda aplicáveis a este 
caso. 

Estes estudos teóricos, no seu conjunto, vêm 
justificar parcialmente os critérios anteriormente 
adoptados no projecto de tensómetros em meios 
exibindo deformações por fluência e que cons- 
tituem o que se designou por «método das inclu- 
sões» [3]. Por outro lado, contribuem para o 
esclarecimento de alguns pontos duvidosos deste 
método, postos no domínio teórico. Assim, pode 
dizer-se que, teoricamente, se comprovam as 
seguintes regras: 


a) Se definirmos rigidez através da relação 
entre força aplicada e deformabilidade (variação 
total do comprimento na direcção de tensão a 
medir) da inclusão, os tensómetros serão tanto 
mais precisos quanto maior a sua rigidez, po- 
dendo esta ser conseguida tanto atraves de con- 
dições geométricas (relação diâmetro-altura no 
caso dos tensómetros cilíndricos) como através 
de condições mecânicas (relação entre o módulo 
de elasticidade da inclusão e do meio envolvente), 
por forma a que se obtenha um nível máximo 
de erro, teoricamente previsível a partir das carac- 
terísticas reológicas do meio; deve ser salientado 
que a rigidez mecânica não tem, do ponto de 
vista teórico, qualquer inconveniente. 

b) os tensómetros deverão ser calibrados no 
interior de amostras do meio envolvente. 

c) será sempre em princípio necessária a utili- 
zação de dispositivos correctores das tensões se- 
cundárias provenientes das variações autogénias 
do volume do meio. 
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A aplicação na prática das teorias recentes ao 
betão necessita, mo entanto, de comprovação 
experimental prévia. Há que considerar o facto 
de o comportamento do material não caber rigo- 
rosamente nos esquemas analíticos adoptados na 
teorização. Um outro ponto a esclarecer é a via- 
bilidade de utilização de inclusões bastante rígidas, 
pois algumas vezes se supõe que elas podem in- 
troduzir concentrações de tensões elevadas, e daí 
produzir plastificações do material, o que iria anu- 
lar a validade da esquematização teórica. Um ou- 
tro aspecto muito importante é o da heterogenei- 
dade do material, constituído por elementos gran- 
des de agregado, imersos numa matriz de arga- 
massa [12] [13] [16] [18] [19]. Considera-se, ge- 
ralmente, que o efeito da heterogeneidade na 
imediata vizinhança do tensómetro pode ser anu- 
lado se este tiver dimensões relativamente grandes 
em face da dos maiores agregados. No caso do be- 
tão em grandes massas, este condicionamento 
obrigará normalmente a crivar o betão envolvente 
do tensómetro, o que cria um novo problema, pois 
as características do betão crivado não serão idên- 
ticas às do betão integral. Neste caso, a zona de 
betão crivado funciona como uma inclusão vis- 
coelástica no seio de um meio viscoelástico. Pro- 
blemas deste tipo são também resolúveis, como- 
já atrás foi citado [16]. Se, portanto, forem apro- 
ximadamente conhecidas as características do 
betão crivado e do betão integral, é possível 
prever um erro máximo, proveniente da criva- 
gem, na determinação das tensões. 

Deve dizer-se que as teorias expostas se apre- 
sentam bastante mais aptas para uma interpreta- 
ção da realidade do que as tradicionais, em que as 
propriedades de fluência do meio não eram carac- 
terizadas. Assim, por exemplo, no caso de ten- 
sómetros alongados, cuja utilização parece van- 
tajosa para a determinação de tensões em ele- 
mentos estruturais finos, ensaios de fluência têm 
por vezes, mostrado uma inesperada progressão 
no tempo das tensões acusadas pelos medidores. 

É frequente atribuir este facto ao efeito da 
concentração de tensões ma zona de contacto 
entre a inclusão e o meio, as quais provocariam 
neste plastificações locais. A existência de tal 
efeito não podia, no entanto, ser cabalmente 
discutida em face de insuficiência das bases teó- 
ricas disponíveis. As expressões agora apresen- 
tadas neste estudo, que fazem intervir as carac- 
terísticas do meio viscoelástico, Em, Bo, Boo, 
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permitem, no entanto, uma separação de dife- 
rentes factores envolvidos nos ensaios, e, por- 
tanto, podem constituir um passo para o julga- 
mento definitivo do problema posto. 

Estes estudos deveriam, contudo, ser precedi- 
dos do esclarecimento de algumas questões pos- 
tas ainda em domínio elástico, como é o caso do 
conhecimento do coeficiente experimental « para 
as diversas relações geométricas de tensómetros. 
Apesar destas dificuldades, as expressões agora 
deduzidas foram já utilizadas [16] para a inter- 
pretação de ensaios conduzidos há alguns anos 
no LNEC, em que se utilizaram tensómetros alon- 
gados em um betão carregado com muito pouca 
idade. A interpretação foi dificultada, entre ou- 
tros motivos, pelo facto de não terem sido esta- 
belecidas as leis de fluência do betão para várias 
idades de carga. Apesar disso, os limitados re- 
sultados que se obtiveram são encorajantes do 
ponto de vista do prosseguimento de estudos 
experimentais acompanhados de uma interpreta- 
ção analítica. 

Para além dos problemas a esclarecer por 
experiências em laboratório, não deve ser esque- 
cido que é extremamente importante uma per- 
feita colocação dos aparelhos para se obterem 
em obra resultados satisfatórios. As dificuldades 
são mais fortemente sentidas no caso de células 
tipo folheto e justificar-se-iam pormenorizados 
cuidados, que nem sempre é possível impor em 
obra. No entanto, poder-se-á talvez ladear a difi- 
culdade fazendo a colocação cuidada em labora- 
tório no interior de pequenas peças de betão 
idêntico ao da obra, as quais seriam postas no 
local pretendido no momento da betonagem. 

São pois, numerosos e de diversa índole os 
problemas a resolver. No entanto, a importância 
da determinação das tensões justifica todos os 
esforços para lhes dar solução. 


APÊNDICE 


Regras Operatórias da Transformação 
de Carson 


Seja f (t) uma função definida no intervalo 
(0,º>) e tal que o integral 


L [ o | ama Fi f Ce-pt ft) dt 


” O 
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exista, sendo o parâmetro p um número, cuja 
parte real é positiva, e escolhido de tal modo 
que o integral seja convergente. 


A função f£ (p) é a transformada de Laplace de 
função f (t). 
A transformada de Carson é dada por 


Pm=pt()=p| ertt(ydt 


À função f (t) podemos chamar original e à trans- 
formada (neste caso, de Carson) f* (p) podemos 
chamar imagem e a correspondência entre as duas 
funções pode simbolizar-se pela notação. 


E (oo f* (p) 


Na tabela anexa apresentam-se as transforma- 
das de Carson e de Laplace de várias funções. 

As transformações de Carson e Laplace têm 
algumas regras operatórias, simples de deduzir. 

Assim, no caso da transformação de Carson, 
temos. 


— Adição : 
tlO+g(lt)o fe ip)+gr(p) 
com o corolário 


K ft)D K E (p) 


— Semelhança : 


E (at) o E ( P 
selo 


— Translacção : 
f(t—a)ne-*Pf*(p) 


— Amortecimento : 


t(e-tp P 
2 = 


f* (p + a) 


a 


— Transformação das derivadas de uma função : 


d" E 
-—— » p"f*(p)—p" E (ot)-p"r-! E (ot) — 
red Mid ae, J—p (0+) 
dr = f 
e. cx qa a - (0%) 


d'f 
em que ão! (0*) representa o limite da função 
p 
: quando t tende para O por valores positivos. 
P 
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— Produto de uma função por uma potência de t: 


n/ E (p) 
d par 
d p” 


ef(t)o (—1)-"p 
— Integração de uma função: 


“t 
| f (7) dz O J f* (p) 
Jo P 
— Convolução : 


E f()g(t—rdr > co (p) g* (p) 


ou 


d t 
Eai f (x É =:t) dtro PP *+ 
E IN (5) e( (p) 8* (p) 


trabalho. 


Para uma exposição aprofundada do cálculo 
simbólico, consulte-se a literatura da especiali- 


dade [27]. 


TABELA 


Este conjunto de regras operatórias permite 
deduzir, com relativa facilidade, a maior parte 
das transformadas de Carson que poderão inte- 
ressar na prática, a partir de tabelas já organi- 
zadas para as funções mais comuns, ou, inver- 
samente, determinar as funções originais a partir 
das transformadas. Na tabela anexa são dadas as 
correspondências mais utilizadas no presente 


IMAGENS E ORIGINAIS DE TRANSFORMAÇÕES DE CARSON E DE LAPLACE 


(p-a) (p-b) (p-a)(p-b) 


p? p 
p?.2Dwp+w? p? + 2Dwp+ w? 


p? 


RE: A | | E 1 
tp-a)(p-b)(p-c) (p-a)(p-b)(p-c) 


p ] 
(p-a)tp-b)(p-c) | (p-a)(p-b)(p-c) 


senh (at) 


eat < ebt 
a-b 


cos wt 


1 Ow sen (wt 1-D2er, ) 
V1-D? 


(1-wt)Jerwt 


dA ePW senh(wtv DijT, ) 


= 


wV1-D2 
sendo T, =arc tg a 


Ww 
« T;zarc tgh 


a(b-cje?!.b(c-ajebt, c(la-bje*! 
(a-b)(b-c)(a-c) 


(b-c)eate(c-ajett «(a-b)eft 
(a-b)(b-c)(a-c) 


— sen (wt) 
w 


e Dwt sen (wtV 1-D?) 


— 1 ePwt senh(wtVD2-1) 


para D:-0 


para 0%) 


para p2. ] 


para D2>1 


/0º.1 


-Dw 


para D=0 


para DX] 
2 
para D'>] 


para D=1 
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AEG-LUSITANA DE ELECTRICIDADE 


M- 4 a 
ta . 
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nai al Ted di 0 Hp 
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, 

- 
em É 
NTE 
2 4 


a imo 


Já na embalagem o nosso transformador 


de intensidade “anão “ 


mostra claramente a sua potência 


E isto, porque «ele» é visivel antes 
que V. o retire da embalagem ou 
que o compre. Sabe porquê ? É que 
a sua embalagem é em poliestireno. 
Basta uma olhadela: Lá estão o 
barramento primário, a chapa de 
características, os terminais e a 
peça de fixação. Poderá imaginar 
outras vantagens : Quando o nosso 
«anão» estiver no seu armazém, não 
necessita de «lhe» remover a sua 
embalagem protectora, à prova de 
poeiras, para saber que é «êle» que 
ali está. O R 02 estará tão limpo 
cómo no dia em que o comprou. 
Nós, AÉG, fornecemo-slo» para cor- 


rentes nominais que vão de 50 A a 
600 A. São múltiplas as suas aplica- 
ções. Também existem para intensi- 
dades inferiores a 50A. A sua 
montagem é fácil, porque o seu 
volume também é reduzido: 0,123 
dm). 


Transformadores 
de medida 
AEG 
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PRÉ - ESFORÇADOS CERÂMICOS 


DATIAL 


patente stahiton 


Companhia Portuguesa de Trelilari 


S.A. R,L. 


ARMADURAS PRE-FABRICADAS MALHASOL 


E 
Tensão de segurança 3.000 Kgf / cm? 


AÇO BI 


/RTON 
NB 


Tensão de segurança 4, 000 Kgf / cm? 
AÇOS DE ALTA RESISTÊNCIA PARA BETÃO ARMADO 


que pela eficiência técnica e económica s que conduzem, são 
um indispensável clemento na moderna construção ctvil 


Para qualquer obra em betão armado cossulte: 


DIAL — GABINETE TÉCNICO 


R. DO INSTITUTO INDUSTRIAL, 18-1,º Dt,º 
TELEFONES 671224/5 e 673785 - LISBOA 


PAVIMENTOS, TECTOS, FORROS, COBERTURAS, ASNAS, SILOS | 


CERVAL 


CERÂMICA VALE.LOBOS.I! 


R. JOAQUIM A. AGUIAR, 41, 5.º-D.- LISBOA - TEL.53 1138/9/0 
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TECNICA XX 


C. D. O. 691.55 


CARACTERÍSTICAS DE ALGUNS GESSOS NACIONAIS 
UTILIZADOS NA CONSTRUÇÃO CIVIL 


RESUMO 


É preocupação antiga do Laboratório Nacional de En- 
genharia Civil o estudo das características do gesso de pro- 
dução nacional e o estabelecimento de especificações de 
qualidade pura este material, de indiscutível interesse para 
o acabamento das construções. 

O conhecimento de frequentes insucessos no comporta- 
mento de estuques em edifícios de relativa importância, 
determinou o início do estudo deste complexo preblema- 
Nesta primeira parte do estudo procura-se definir o nivel 
de qualidade des gessos produzidos na metrópole em com- 
paração com o exigido pelas principais normas estrangei- 
ras e salienta-se a necessidade de ensaios sistemáticos que 
permitam aquilatar da constância da sua qualidade e con- 
tribuam para o estabelecimento dos limites a prescrever para 
as suas caracteristicas quimicas e físicas. 


1 — INTRODUÇÃO 


O gesso é, como se sabe, um ligante hidráu- 
lico obtido pela cozedura do minério também 
designado por gesso ou gesso bruto, essencial- 
mente constituído por sulfato de cálcio hidratado 
com duas moléculas de água — Ca SO. 2 H;0O. 

Na sua forma mais corrente, o gesso cozido 
apresenta-se como um hemihidrato, quer dizer, 
contém meia molécula de água por molécula de 
sulfato de cálcio — Ca SO;. */» H:O. Moiído e 
amassado com água faz presa, reconstituindo-se 
o hidrato com duas moléculas de água. 

O gesso é conhecido como material de cons- 
trução desde remota antiguidade e já utilizado 
muito antes da cal nos países de clima seco, como 
a Abissínia, a Caldeia, a Fenícia e o Egipto. 

A fraca resistência mecânica e a sensibili- 
dade à humidade fazem com que a sua utiliza- 
ção na construção civil, entre nós, se limite quase 
exclusivamente ao acabamento de paredes e tec- 
tos e ao fabrico de placas para a execução de 
divisórias. 

Como material de acabamento que é, a sua 
qualidade tem influência relevante não só 
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SYNOPSIS 


Itis an old intention of the Laboratório Nacional de 
Engenharia Civil to study the characteristics of qypsum 
plaster made in Portugal and the preparation of quality 
standards for this material, vf doubtless interest for buil- 
ding finishes. 

This complex problem began now to be studied due to 
the unsatisfactory behaviour of plasterwork, frequently 
observed in relatively important buildings. In the first 
part of this study the quality of qupsum plasters manufac- 
tured in Continental Portugal as investigated in compa- 
rison with the qualimy requirements of foreign standards, 
and attention is drawn to the need of systematic testa for 
Judging if the quality of the products is uniform and for 
giving data which can be used for setting down limit 
values of their chemical and physical characteristics. 


no aspecto das superfícies e na sua conservação, 
mas também na durabilidade e efeito decorativo 
das pinturas que sobre elas se apliquem. 

O comportamento dos revestimentos estucados 
depende dum conjunto de factores entre os quais 
sobressaiem : além do próprio comportamento das 
bases, a qualidade dos materiais empregados, as 
técnicas de execução dos trabalhos e as condi- 
ções do ambiente em que decorrem o endureci- 
mento e a secagem desses revestimentos. Por 
estes motivos, é preocupação antiga do Labora- 
tório Nacional de Engenharia Civil o estudo da 
normalização das características destes materiais 
e das técnicas da sua utilização. 

O conhecimento de frequentes insucessos no 
comportamento de estuques, em construções de 
relativa importância, determinou o início deste 
complexo estudo que visa, naturalmente, a obten- 
ção do conhecimento tão completo quanto pos- 
sível dos gessos disponíveis no mercado nacional. 


2 — AMOSTRAS DE GESSO ESTUDADAS 


Durante os trabalhos preliminares do planea- 
mento do estudo, supunha-se haver, na metró- 
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pole, cerca de 10 fábricas de gesso, com capaci- 
dades de produção muito diferentes. 

Algumas das amostras que se consideraram 
neste estudo foram obtidas em duas procuras fei- 
tas no mercado de Lisboa, em meados de 1962 
e de 1964, e outras foram remetidas ao Labora- 
tório por uma fábrica algarvia que tinha reque- 
rido a verificação da qualidade da sua produção 
de gesso de acordo com as exigências do decreto 
n.º 42 941, de 23 de Abril de 1960. 

Da primeira procura feita no mercado de 
Lisboa (!) em 1962, resultou a colheita de amos- 
tras de gesso produzidas por 4 fábricas diferen- 
tes, enquanto na segunda, se obtiveram apenas 
amostras de duas fábricas diferentes. Tal facto 
resultou das quatro primeiras fábricas se terem 
agrupado e passarem a embalar todos os seus 
produtos com um único nome comercial e de se 
ter então encontrado gesso de outro produtor. 


QUADRO I 


Identificação das amostras utilizadas 
neste estudo 


Amostra 


Fábrica Designação 
aa — ou do 
N.º Data da colheita marca Lesso 
1 A | Janeiro 1962 ii | ia 
1 B | Fever. 1962 | 
2 B Extra 
3 C Estuque 
Junho 1962 
4 D Paris 
õ E | Francês 
6 Outubro 1962 | 
A 1,4 
7 Janeiro 1963 
8 1. 
ein F 
9 9 a 
—— | Abril 1964 ia 
10 | 4 a 
na G 
1 2 À 


(!) Houve a preocupação de colher o produto directa- 
mente no mercado e não de o solicitar às fábricas produ- 
toras, 
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Quanto às designações das várias qualidades de 
gesso existentes no mercado, nota-se estarem eiva- 
das de uma falta de uniformidade que seria conve- 
niente eliminar. Assim, as várias fábricas usam 
as seguintes designações indefinidas, para identifi- 
carem os seus produtos: 

— Extra 

— Estuque 

— Paris 

— Francês 

— 1.2 qualidade 
— 2.2 qualidade 

No Quadro I identifica-se cada amostra estu- 
dada e, por razões facilmente compreensíveis, 
omite-se a sua proveniência, mas mantém-se 
a designação com que foi obtida no mercado, a 
fim de verificar a existência de qualquer correla- 
ção entre estas designações e as características dos 
gessos respectivos. 


3— CARACTERÍSTICAS EXIGÍVEIS PARA O 
GESSO COZIDO 

Não estão ainda especificadas entre nós as carac- 
terísticas de qualidade exigíveis para os diversos ti- 
pos de gesso cozido a fim de que se comportem sa- 
tisfatóriamente nas utilizações a que se destinam. 

Por isso se recorreu à comparação das caracte- 
risticas dos nossos produtos com as exigidas 
pelas normas estrangeiras abaixo mencionadas, 
o que implicou a utilização das técnicas de exe- 
cução dos ensaios indicadas por essas normas. 
Por esta mesma razão não se deu preferência aos 
métodos de ensaio estabelecidos pelas normas por- 
tuguesas disponíveis, em fase de inquérito nessa 
altura. 

No quadro Il indicam-se as normas portuguesas 
referidas bem como as que foram publicadas pos- 
teriormente como definitivas. 

Para facilitar a comparação da qualidade dos 
produtos nacionais em estudo com a exigida 
pelas várias normas estrangeiras consultadas, apre- 
sentam-se a seguir os resumos destas normas no 
que se refere sobretudo a tipos de gesso idênti- 
cos aos ensaiados. 


3.1 — Norma Francesa 
NF B12-001 — Aoút 1946 «Plátres de la Region 
Parisienne» (*). 


(7) Esta norma foi substituída posteriormente pelas nor- 
mas NF B12-300, NF B12-301 e NF B12-401; a este as- 
sunto refere-se também o fascículo de documentação FD-B 
N.º 12-402. 
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QUADRO II 


Normas portuguesas sobre gessos 


Normas existentes antes da execução dos ensaios 


Normas actualmente existentes 


| | Designação ai Designação 
Identifl- para (1). pesos Data (1) 
een | | jeral | Específica " | Geral É Especif ica 
| nm e 
I-313 | 9/60 | Gessos | Terminologia | 
i ssa mu a — a 
| Características, 
| transporte, 
I-314 | 9/60 Gesso cru identificação e | 
| | recepção | sm 
MRI FND PO 1 o 
| | | 
1-315 | 9/60 | Gessos Colheita de amostras | NP - 315 29/11/63 |s o | eiiainaigida 
| RR Eno o 
E | + 
| | Capaci idade de absor- NP -316| 29/11/ | Classificação, carac- 
A ' - 11/63 |! Gesso bruto | 
1 - 316 | 9/60 | Gesso cozido] ção de água | [a e recepção 
| 
| | 
I-317 | 11/60 | Gessos | Humidade NP - SL | era /63 | Gessos | Colheita de amostras 
ms | ssa | : 
| | Determinação da 
| SS0s Água combinada NP -318| 29/11/63 Gesso 
I-318 | 11/60 Gessos |, E / | água de presa 
o 
—— | E " Ki 
| à Princípio de presa | Determinação da hu- 
' | o 19 | 29/11/63 Gessos 
+ PAM a1/ so pfrso on e tempo de presa E á | midade 
e q a a tm: = a 2 
1-320 | 12/60 | Gessos Teor em anidrido | np. 390, 29/11/63 Gessos | Determinação do teor 
carbónico | em água combinada 
| | Teores em resíduo in- | | Net . 
us | solúvel em óxidos de am Gi ; u eeião pen cmo qo 
I-321 | 12/60 Gessos ferro e alumínio, de -321, 4/1/64 ezso pes (8) o presa e 
cálcio e de magnésio nl a 
E E | 
| oca. Determinação do teor 
pe ] Teor em anidrido 20/11/63 Gisiós anidrid REA 
I - 322 | 12/60 Gessos ESTA NP - 322 | 29/11/ E=s0s em ani ú o carbó- 
| | | nico 
RE | | 
| | Determinação dos teo- 
Teor em cloreto | Padligi ceia 
I-323 | 12/60 Gessos de uôdio NP-323 29/11/63 Gessos lúvel, em ferro e alu- 
| | | | mínio, em cálcio e em 
| | | magnésio 
i | 
L-324 | 12/60 E Apresentação dos re- - 324 | 29/11/63 | diataica Determinação do teor 
| sultados analíticos | | em sulfato 
ri | 
| NP -325 29/11/68 isaanca Determinação do teor 
| em cloreto de sódio 
pç ni E | 
P-326 14/8/62 Gessos Apresentação dos re-| vp. 996 4/1/64 Gósioá Apresentação dos re- 


E - 369 | 


| Gesso cozido 


sultados analíticos 


uso e e 


Granulometria 


sultados analíticos 


ca) Pol essreseem À La! 


| | NP - 379 | 19/12/64 | Gesso 


Granulometria 


(1) As datas indicadas são as de publicação dos inquéritos (1) no Boletim de Normalização; as de homologação 
das normas provisórias (P) e as de aprovação, por portaria, das normas definitivas (NP). 
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